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ABSTRACT 



An optimal method of filtering radar data is analyzed and experimented. 

First the theoretical bases of Kalman's approach to filtering 
problems is presented. Then the suitability of Its conclusions to 
the filtering of radar data is studied. The form for the filter needed 
in such a case is derived and, finally, a simulation is made attempting 
to show the real usefulness of Kalman's filter. 

It is concluded that the method considered presents a very prom- 
ising set of features which may lead to advantageous practical use. 

The filtering process Is based on the recurrent solution of matrix 
formulas. To accomplish a reasonably fast computation of the filtered 
values the solution was programmed in the FORTRAN language and run in 
a 1604 CDC computer. 
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I. Introduction 



We will be concerned in this work with the filtering of radar data. 

To our knowledge the solution of this specific problem, in a major- 
ity of cases at least, was never attempted with the help of a deep theo- 
retical foundation. So far tne "«> Tracker" seems to have conveyed 
almost all the possibilities of practical processing of radar data. 

With the new approach to the general problem of filtering announced 
by Kalman [l ] a new frame of practical processing stood up to the engineer. 

The " «<,|3 Tracker" has no inherent means or need of utilizing the 
known statistics of the targets to be tracked in any mathematical or even 
rational way. 

The Kalman approach starts out by quantifying the known parameters 
of the signal generating process. Broadly speaking, the class of prob- 
lems attacked by Kalman consists of filtering and predicting the output 
of Gauss-Markov random processes when its parameters and statistics are 
known. The method described by Kalman being particularly suited to 
machine computation was immediately explored in more than one field 
of possible applications. 

One branch of the work done with the Kalman filter occupied control 
engineers. Another one gave results in the tracking of artificial 
satel I i tes . 

Rauch [2] and Magi 1 1 followed Kalman with the analysis of 
cases in which the process parameters are random variables. 

It will be shown that Mag ill's Technical Report is completely 
suitable to provide a theoretical basis upon which a radar data filter 
can be built with fairly good results. 



2. Statement of the problem. 



The processing of data provided by the. various typos of radars 
can be readily implemented with the digital computer. 

The processing device must be able to allocate the echoes received 
to one of the various targets being tracked, to decide on the appearance 
of a new target, to give numerical significance to the distance, azimuth 
and other characteristics detected by the radar, to store sequentially 
the information being obtained and to provide all the data stored to 
any unit or entity allowed to obtain it. 

This data will be considered to be present and available for 
filtering purposes. 

The data being stored needs to be handled in such a way that the 
various tracks maintained are filtered as accurately as possible. It 
is also considered that a prediction will have to be made and trans- 
mitted to some other system. 

The characteristics of the data provided by the computer storage 
for an eventual solution of the filtering problem are, obviously, of 
great importance. 

Admitting that the allocation of echoes to targets was done 
previously with total confidence, we do not need to consider but the 
filtering of one target track. 

This way for each target we will have stored data with the 
following characteristics, 

a) Numerical values assigned to distance, azimuth and altitude. 

b) Times at which the position was detected. 

c) If relevant, data on the position, speed and course of the 

vehicle where the radar is mounted. 
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Adding to this data collected in the machine storage all the char* 
acteristlcs of the target defining its dynamic model are also necessary 
to improve the tracking efficiency. We may specify at this point that 
it is absolutely necessary to have some "a priori" degree of knowledge 
about the dynamics of possible targets when the first point of the track 
is treated in the filtering process. As will be seen, the filter will 
not start its recurrent computation unless previous information on the 
characteristics of the possible targets is available. 

We need then: 

d) A dynamic model capable of describing the motion of the tracked 
object. 

e) Some statistics of the inputs to that model. 

Finally, the source of the noise to be filtered out has to be 
considered. Noise will be introduced by the radar itself, by the 
processing prior to the filter action, and by any anomalies in the 
propagation of the electromagnetic waves. 

Let us consider then: 

f) Statistics of the error introduced by the equipment. 

With this data available the filtering process will atm at 

reproducing the real track of the target, eliminating errors added 
during the acquisition of data. 

We will see that this objective is not completely achieved but 
that the mean square error of the positions will be reduced to a 
minimum, with the help of the Kalman filter. 

In d) we referred to the need of a dynamic model capable of describ- 
ing the motion of the object being tracked. More precisely we can say 
that this model is going to include or use, In one way or another, all 
the data gathered by us or available in storage. We may even say that 
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this dynamic model is central to the solution we are seeking. With 
this in view let us proceed to formalize the model and obtain a more 
rigorous concept of the problem presented to us. 

The signal or message and the noise precesses are given by, 

5 ( t-h/) - $ s (t)$(t) + 4 ( t+t) 4 ( 

yd) * r> sfrj 

and, 

W(t+>)~ Ct) u(t) + A (t+j) u w (t+/) 

n (t) ~ h r (t) w (t) 9 

Where the various symbols have the following meanings: 

SCtf is the present state vector of the message process, i.e., 
a column matrix formed by the target components of distance, velocity, 
acceleration, etc., on some coordinate axis (which we assume to be the 
X-axis from now on). In chapter three we will indicate how the other 
components can be treated by the same filter. 

is the transition matrix of the message process, interpreted 
as the matrix by which we must multipy the present state to obtain the 
next, as will clearly be seen in the end of this chapter. This includes 
the knowledge of d) . 

U,(t) Is a vector that includes all the possible inputs to the 
message process, as wind gusts (which we will consider as being able 
to change instantly the position of the target), air speed and varia- 
tions of acceleration as caused by the pilot's action. As later will 
be seen, these inputs have to be white gauss i an random variables in 
order that the theory can be applied. 

J) 9 (tj is the distribution matrix of the message process and 
determines in which way the inputs affect the various state vector 
components . 
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y(t) Is the message output which may consist of one or more of 
the state vector components, 

Is the output vector of the message process which determines 
what state vector components are going to be used in the model output. 

(frj Is the transition matrix of the noise process. 

is the white noise process vector. Includes the knowledge 

of f). 

ti(fc) is the noise process output. 
hfi) Is the output vector of the noise process. 

The central problem of this investigation can now be clearly defined. 
The information received by the radar and sequentially stored is the sum, 

z(t]~y(t)+n(tf , > 

where y(t) represents the true distance and r>(t) tli i , .ir added to 
that. We will seek a procedure of machine computation to obtain y(t] 
at each instant t- 1,2,3,.... and, in a more general fashion, to obtain 
some y(i~)j-€y ( where the error e /(V is minimized in accordance 

with the mean-square error criterion. 

* 

Notice here that as we are treating the X-axis components this 
will have to be interpreted as distance along the X-axis. This apparent 
confusion could be avoided but we decided to maintain completely Magill's 
notation to provide some unity between his report and this application, 



3. Consideration of the theory 



Later on we will study Kalman filter's adequacy to radar data filter* 
ing and will concentrate in the areas of application. But, obviously, 
a clear understanding of the theoretical results to be used has to be 
provided. This Is the present chapter's objective. 

The elements required for a solution may be concentrated in the 
four equations presented in the last chapter which describe the message 
and noise processes. 

A simplification of this can be made and consists of joining to* 
gether the message and noise processes, reducing these to a single 
process which has as output the sum ~z(t) = '/&+ n (&) This is accom- 
plished by augmenting the state vector. 

According to this we now define: 




vector of the observable process. 




transition matrix of the observable process. 






output vector of the observable process. 



Defining also, 




covariance matrix of the inputs. 



!>(*)& 
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The model for the process which generates the observed signal 
may now be reduced to the set of equations: 

X (&+/)= <t>(t) x(t) t 2>(tn) U(t+t) 

2 (&) = > 
where U(t+7) 

Even if this group of definitions casts some shadow on the physical 
meaning of the parameters we will not try to clarify it right now. The 
next chapter will provide full understanding by means of an example. 

We may now represent the whole process in a block-diagram form. 

Fig. 1, where the parallel lines represent flow of vectors instead of 
a scalar value as in the line at right for 2-(&J. 




1 > 



Fig. 1 Model of the observable process 

- 4 , 



*As Q(tJ is a diagonal matrix, P is equivalent to apply 



the exponent — - to all non-zero elements of this matrix. 
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Let us proceed now to look at the solution of the filtering problem. 

The objective is to find a recurrence formula to obtain for each 
sampling point the best estimate of the state vector, given both the 
previous one and observed state vectors at that sampling point. 

The best estimate of the state vector we need is by definition, 

KCt/t-O = £ jx(t) /Zti) , 

Where represents all the available data up to and including time t - /. 

Z(t/N)i Z(tJ-2(t/t->) , 

where z(t) is the observed signal at time t and is the best 

estimate of the observed signal at time t given the previous best 
estimate. It is possible to derive [3j a recurrence formula, 

X (ir/t) - * ( t/t-f) + t/t-t) , r , 2, . ... , (/) 

where K(t,t) represents a variable gain vector. 

We now observe that the filter produces the best estimate at time t 
given the data collected up to and including that time by means of the 
following operations, 

a) Evaluates a best estimate at time t given data collected up to 

and including £•/ . This is simply done by multiplying the transition 
matrix j by the best estimate • 

b) Corrects X ( t/t-l) by the addition to It of a term containing 

the information received at t . This term is the product of a weighting 
coefficient - the gain - and the error in estimating the observed 

signal given data collected up to and including . 

From what has been said we understand that the filter is essen- 
tially formed with a model of the message process and the added feature 
of being able to correct the updating of its output with the information 
received. 
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The recurrence formula (1) is the solution to the filtering problem. 
Of course the means of updating the gain matrix is still not uncovered. 

The derivation of a recurrence formula for this matrix starts out 
considering that the best estimate can be obtained by minimizing the 
trace of the following matrix, 

P(t/t)i £ {x (t/tj y T (t/t) J 

where X (tit) ” X(t) -X (t/t) , 

X (t) being the true state vector of the message and noise gener- 
ated; 

The matrix is the covariance matrix of the estimation error 

and, obviously, gives an indication of the error value at any t' . 

Further defining, 

r * (£/£-') = P (tit-/) tv (&)- M''jz(t/t-fj ( 2 ) 

it is possible to obtain 

M(t, t) - P(t/t-/)m(t) (3) 

(T 2 (f/t-l) 

A recurrence formula for /) must still be obtained. Again 

Magi 11 derives it based on the assumptions previously made and reaches 
the result, 

P(t+l/t)~ 0(t)[l-k(t, t) m YtJ ]P(t/t-t) 0 T (th W* Yt) fid 

The recurrent procedure is now completely obtainable, provided the 
Initial covariance matrix P(t u / previous knowledge,) is given. This is 
done by sequentially evaluating ^(t/tJ,P(t¥-//t) , etc.. 

With the equations we have introduced, a filter block-diagram is 
possible. 
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Fig. 2 Model of the optimal filter 
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k. Applications 



In this chapter we intend to analyze more closely the assumptions 
made before getting the results of the last chapter and test them in 
the problem of filtering radar data. We will then decide on the 
applicability to our problem of the theoretical results presented before. 

As It was said, all Hagill's derivations are based on the fact that 
the processes involved are either Gauss-Markov or can be obtained from 
these by application of a linear matrix multiplier. Can the state 
vector of any flying target, whatever dynamic properties it may have, 
be considered a Gauss-Markov process? In fact, given any instantaneous 
state vector, it is always possible to consider some other previous 
time when the definition applies, i.e., 

f=>[x(t+0/x(t), t[xf t+0/* > 

where K is finite. 

In a physical sense we may say that if we want, at any instant t » 
to obtain the state vector at t+i , we may discard the knowledge we have 
of the target's state vector before the instant t-/C . For example, If 
we consider an aircraft at t , clearly it is seen that the knowledge 
of all its positions since the flight started is normally irrelevant to 
determine the next one but we are almost certainly interested in the 
last two or three. Of course the number fC will vary with the type of 
target considered but in any case is a finite number. 

In Chapter two we said in d) that we would need a dynamic model 
for the tracked target. This model and the information of e) (some 
statistics of the inputs) and f) (statistics of the radar error) were 
all considered in the formulas then presented. 



11 



Let’s analyze first the message process , 

s(tn) =■ (p s (t) s (t)+ 4 

y(ty • r T (t)s(t) 

If we consider In the state vector only distance, speed and 
accel leratlon and this way eliminate the possibility of a third deri- 
vative In the motion we have. 











%(t) 




■£,<W 




%(&h) 


ac 


pjt) tjt) tjo 




%(*) 


+ 


4 


W 


to to 

















Simplifying the notation and multiplying matrices we obtain. 



* A s , * <P a S ‘ * K s * * 4» u s, * 4,2 4, * 4, 

The first three terms of the second member of the equation repre- 
sent the formula to obtain the distance at some Instant t*-i , the dis- 
tance, velocity and acceleration at t being known. The components^, 
<p n and could be easily gotten through the Taylor series expansion 
and as we know they are 

#,=' > r , 4 , , 

where 7~ Is obviously the sampling period. 

The other components of the $ matrix are easily obtained In the 
same fashion. 

We can now give full significance to D s O s if we recall 

the definition of the U$(t+i) matrix given before. It is simply the 
way of correcting the next value of the target’s state vector by adding 
the possible random inputs. For example, we may obtain the target’s 
next velocity value if we know the present one, the acceleration and 
the sampling period but notice that we still have to consider the intro- 
duction, during the sampling period, of some wind velocity and then we 
have to complete the value of S z (&)+ TS s (t) by the addition 
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Of 44 (t+l) . 

The random Inputs Uyft+O and U$ 3 (t+/) account for Instantaneous 
changes of position and for variations of acceleration Introduced by 
differences In the engine's thrust. 



All put together in the matrix equation we have, 



S, (M) 




' / r ±T Z ~ 




s,(t) 




too 




^ 5 / &+0 


S, ft*/) 




0 / r 




S 2 (t) 




O/O 




U u (tn) 


S /*,) 

_ - 




0 0 t 




1 




0 0 / 


1 


y a (** 0 m 



One feature of the random inputs when considered together is that 
one of them may substitute the entire vector. Consider, for example, 

• Through this acceleration Input we can account for any 
variation in position or velocity which instantaneously may happen. 

Looking now at the second equation of the message process and 
considering y(t) as a vector containing only distance, 



L'Ct)] = 



I 

0 

o 



In the noise process we have a different case. We consider the 
noise at any sampling time to be independent from any other noise value. 
This makes all the elements in null. 

Let us make it clear also that the state vector has only 

one component which is, obviously, the error introduced In the measure- 
ment of distance, the only one made. The Input U»(&h) is then the 
actual error \/tt+/) . 

We have then, 

[w(tH)~]= [°J[ W(tj] +[l] [(S» (t+!)] j (*) 

and, 

[»(*)] - [0 [*(*>] to 
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A point we still did not look at is the gaussian quality of the 
random inputs. Can we admit that the acceleration imparted to an 
aircraft when considered in discrete amounts constitutes a random 
variable with gaussian distribution and zero mean? If we think of 
all the aircraft and all its flights we, naturally, would be ready 
to accept that. And actually, this is what happens once the filtering 
must be appropriate for any aircraft that can appear in the radar scope. 
This has to be kept in mind when values are going to be given to the 
inputs. We will return to this problem in the next chapter. 

Having formed the message and noise processes we may now proceed 
to construct the model for the observable process, as explained in 
Chapter 2. 

So, 



x(£)* 



£,(*) 

Sx(t) 

*3 CO 
wit) 
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and with simplified notation, 



Q Ct) 



oS 


^S/^SZ 


^ i / 




Ue/ ^Si 


^55 






U&! Vf* 




&S3 ^So> 



As the random variables 
see that , 

£ ^ > 

and, 



£ )<J S , 0 K }= Co/fa, U t )*0 



have zero mean, we 



s i nee , etc. are uncorrelated. 



(Jo) '* O O o 

o o o 

O O 0 

° o o c** ru s*y* 

<ro s , ooo 

0 iTti, 0 0 

O O Gil, O 9 

O O o CTi*_ 

where (T represents the standard deviation. 

With the matrix just obtained for &Ct) we are able to use the 

recurrence formulas presented in the proceeding chapter. 

But we still need the initial covariance matrix without which we 

cannot start the computation. 

It can be observed very easily that at , i.e., when the first 

measurement of the observable process is made, the correlations between 
distance and speed errors, distance and acceleration errors, etc., are 
all null. We have then. 






/ 0 
o t 
O 6 
O o 



0 o 
o o 
/ o 
0 / 



//>r ext' 0 us kr>ow/e</?<?) s 



» kr'K, o 0 O 

O &rj( 2 o ^ O 

O o o 

0 0 0 
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Of course the diagonal values must be based on the information 
about the target which we have even before we begin tracking it. To 
avoid a purely conjectural evaluation of the variance we may simply 
discard the three first received positions from the radar, as far as 
filtering is concerned, and make use of them to estimate the velocity 
and acceleration of the target. At £. we will have then the position, 
the velocity and the acceleration of the target with errors the variance 
of which we can determine. 

Let us see what are those variances. Supposing that we know the 
variance of the error in the distance measurement? 



because the error in the distance measurement is the error introduced 
by the radar. 




If we take X/fe+Z) as the observed value, then we have. 



l&r X, (t* +2) «■ i&r 0^" 



We want now )/<*<* X * (to +2) . Using the last two points we have, 



using the same assumptions as above. 



X 2 (t^Z) _ X,fto-h2)- X,(to^/) _ _ 

7 - 7" 

- XiCfc+z) - K.C k+2) _ /*> ( £>+/)-*, C& +/)1 _ 



r 




Let us obtain now 



(dr ftb+2) 



The three points we have permit us to write, 



Xs (to *£) = (t~ 0 +2) — X 3 (t~o i-2) 3 * 

a? - 2 Ql QtkUl _ *3 ~ X t %*■/) a 

T T 4 „ 

y,(&+2)-x,(tiH) x(M-Xfh) $,<L+2)-X t (M Xfc*)-X,fa) 
t ~ t t ~ r _ 



T . r 

— X/ (ib*-z)~2x t (%i-f)^X,fh ) K&t2± Z&L Qktil ±*£ (2° l 



Now 



]/dr X? (toi-z)m -L^ jyjnZ t fe+$+4terX('6 > i-/)+i/'irX t (h)J =. 



= .£ ^ OX 

T* 
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5. Computation of the filtering process 



In the beginning of this work we showed our interest in the 
eventual machine computation of the recurrent process. It is now 
time to look at the procedure introduced in the last chapters and try 
to evaluate how effective the machine computation can be made. 

in accordance with what we mentioned in chapter three the 
formulas below constitute the recurrent process, 

O’Yt/t’O* n> T (tJ P(t/t'/)f77 ft) 
i ( t/t*)* m ( t-/)k (t-i/t-/) 

The program for this computation can be done very easily. We 
have in the next page a very simple flow-chart that illustrates 
this array of computations. 
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Fig. 3 Flow-chart for recurrence formulas 
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The problem of how to start the recurrent process has still to 
be considered. We said before that we could obtain the value of the 
initial covariance matrix by simply taking the three first points of 
the raw track, evaluating with them the position, velocity arsd accel- 
eration at the third point and, at the same time, calculate the initial 
covariance matrix of the state vector thus obtained. 

The formulas then mentioned are repeated here for convenience, 

x, (b+i)=y(t*+2.) 

X* (tb+Z). Vfr 

T 



X, ( &tz)= z H 

7 " 



X46 h+2) « 0 

X/ (jto /■<? J m tfdf* C/st%f 

l/£r* Us*/ 

i/dr X? (to+2) • — iXnr Ch/ 

7 "* 

Xdr U Sw 



Considering the computations above we obtain the flow diagram 
on the next page . 
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Fig. 4 Flow-chart for the optimum filter 
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This way we use the initial covariance matrix at V«+Z (JT^y) 
to obtain the next one at and use the observed positions at to^2 

to obtain a best estimate at . Our filter will then begin filtering 

at and including t^+3 . 

To avoid too long statements which would add only to confusion in 
the program we divided the computation of the covariance matrix in 
parts as it can easily be seen. 

Another point in need of consideration is the way the program is 
implemented without the help of sub-routines, even if some matrix 
operations are repeated several times. We thought that some computing 
time would be saved by eliminating the call of sub-routines. This way 
there is no great compilation time spent and probably the compilation 
done comes out more efficient. 
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6. Example of a simplified simulation 



In this chapter we are finally ready to Implement the techniques 
which have been studied until now. Ideally that should be done in a 
real case. A seemingly feasible way would be to obtain radar data from 
maneuvering aircrafts, project the points in the X,/ and 2? axis and 
try to filter the three tracks thus obtained. This way of Implementing 
our methods Is, however, not too flexible. 

A better way Is to create a track, add some noise and try to 
filter this out afterwards. This approach will provide us with the 
possibility of varying the track at our will and thus permit the 
evaluation of the filter's suitability in various conditions. How 
much the results can be relied upon depends, obviously, on the simu- 
lation's fidelity. This we will try to assess in the next chapter. 

As we have said before, to study our filtering problem it is 
irrelevant to consider more than the projection of the target's track 
on one of the coordinate axis. The procedure for the other axis is 
exactly the same. 

Limiting ourselves to the projection on the X-axis we can have 
three different situations - constant velocity, constant acceleration 
and sinusoidally varying distance. These correspond in three dimensions 
to constant velocity, constant acceleration and circular track. Any 
other track of Irregular form can be reduced to a group of circles with 
accuracy as high as we want. 

It seems then only reasonable that we study tracks where these 

three possibilities happen in different orders and with different 
values. From here we may eventually conclude that the filter provides 
acceptable results in any conceivable simulated situation. 
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Furthermore it seems 1 eg 5 call 5 y acceptable to conclude the following: 

a) The filter will yield good results with any constant velocity 
given ftp that Is the case for the lowest and highest expected velocities 
simulated . 

b) The filter will work In good conditions for any constant 
acceleration given that It works for the highest possible and fulfills 
the requisl te a) . 

c) The filter will definitely accept any simulation and yield good 
results if, furthermore, it gives acceptable results with the tightest 
circular maneuver possible - In the sense that the acceleration on the 
X-axis reaches the highest value possible with the speeds available. 

At this point we will have to decide what types of targets are 
we going to be handling. Let us confine ourselves to airplanes. 

To simulate a track we will suppose an aircraft goingrthrough all 
the situations presented above. 

Furthermore it seems desirable to have all the combinations of 
those pieces of track that seem to introduce different degrees of 
difficulty. 

The tightest maneuver we want to explore Is the one which gives 
the greatest maximum acceleration of the motion's projection on the 
X-axis. Ot is easy to see by Fig. 2-29 in one of the references DO 
that the highest speed and greatest bank angle will give the largest 
value for acceleration. Considering that the maximum possible man- 
euvering speed for the aircraft is 1000 knots and the banking angle 
70 degrees, we have that the projection on the X-axis of a circle 
described in the X~Y plane is, 

X * /€ Cos x 

where ft is the radius of turn in miles, the rate of turn in degrees 
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per second . 

Extrapolate 8 ng 5n the tables for 1000 ml we have, 

R = 32.000 ft -- ft* 6*4 m/ 

/ O/OOO 

/iot-s 3*6 X = • 05237 rad /sec 

/&o 

Then the maximum acceleration ls» 

R'/ rotf>f = R'(frt) 2 *? 6*4-x{oSZ37)* rot J sec 2 - 

«. C>4x* 0627? f7?// zee* *■ -O/Zf m/ /z&C* 

Let us assume that In a straight line the maximum acceleration 

possible 5s 1.805 g"s. This gives, 

f t 305 X32*2 ft/zec 3 - &tlL rr?/Aec * 9.63 X./o*mt Jsec* 

6 k / o 5 ' 

To have the largest range of speeds we wi 1 1 consider that our 
aircraft can fly from 10 to 1000 mi /hr, which is unreal but will give 
greater scope to our results. 

Let us see now what types of tracks will be worthwhile to exper- 
iment with. Below are the different combinations which seem to present 
different degrees of difficulty to our filter. 

a) Lowest speed, acceleration to highest speed, highest speed, 
circle at highest speed. 

b) Lowest speed, acceleration to highest speed, circle at highest 
speed, highest speed. 

c) Onverse sequence of a). 

d) Inverse sequence of b) . 

Another factor that may be important 5s the time during which the 
aircraft will stay at the lowest and highest speeds. 

Two problems are still not solved. One 5s the addition of noise 
to the track we are going to create. The other is the sampling interval 
which corresponds to the radar rotation period. 
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The noise will be considered to have a Gaussian distribution with 
a variance of four miles constant for any distance. Of course this is 
a rather strong simplification and in reality the variance will depend 
very much on the distance. We do not think that it is too important 
to deal with this problem right now. In the next chapter this will 
be discussed again. 

The sampling period will vary between some limits in a real 
situation. But, again, that fact does not seem to be important in 
this simulation. 

The programming of the radar data simulation is very clear and 
can be analyzed in the Appendix. 

First, a time array is generated, Then, for each piece 

of track we proceed as follows: 

For each sampling point the noise is generated with the RNDEV 
sub-routine which provides us with random numbers obeying to a 
gaussian distribution with variance of unity. DEV is the output 
number, which we multiply by n/coe to get a variance of valve COE. 

Then we evaluate the signal output vector using formulas of elementary 
physics. Finally, the distance is added to the noise giving the input 
to our f 1 1 ter. 

A last question Is what values to give to the $ matrix. First 
of all we are going to assume that all the random inputs are null 
except the acceleration. Then we will try, by the use of several values 
for , to obtain one which fits best the whole filtering for all 

the cases we are going to experiment with. 

The values that finally were reached are printed out In the Appendix, 
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To help In the evaluation of the filtering quality a smoothness 



index was defined as, 

smooth - tirfrvr) y 

*»(*rr) /e 
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7. Concl us Joins and acknowledgements 



The (results obtained indicate that the Kalman filter can be used 
in the filtering of radar data with some success. 

The cases taken as examples did Indicate that in a wide range of 
real situations the filter will provide useful filtering, as far as 
the simplifications made are not the reason why the filter worked. 

The constant sampling period which was used could have been varied 
along the track without any difficulty. To prove this we only note that 
for each computation made after the arrival of a new radar position the 
value used for T 5s not related in any way to the last value of 7”. 

The radar error may have a bias depending on the distance. This 
bias has to be experimentally determined and its value subtracted from 
the range measured. If not it will eventually introduce a bias in the 
filtered values. 

The Gaussian quality of the radar error 1$ a point of doubt. Only 
in the presence ©f a definite radar system can that be asserted. How- 
ever in a general way we observe that all components of the radar noise 
seem to have a Gaussian quality. 

The value we should give to the Q matrix 5s critical, ideally, 
for each type of target (helicopters, airplanes, missiles, ships, etc.) 
we should have a Q matrix. The values for can be chosen experi- 

mentally. The filter will work with less error when the variance of 
random inputs 5s small, as it would be the case of a missile when com- 
pared with an aircraft, especially if we consider a ballistic missile. 

The components of the motion ©n the three coordinate axis can 
obviously be evaluated in succession and with a normal speed computer 
the computation will be completed before the next sett of data arrives. 
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Before terminating we would like to express our appreciation for 
the help given during our work by Or. Harold A. Titus and LCDR Frederick 
0. Jardine. The first was always ready to clarify, challenge and 
ultimately give assurance to our ideas and the second worked very 
closely with us during the study of the theoretical foundations of 
the Kalman ft I ter. 
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APPENDIX 



OPTIMAL FILTER PROGRAM 



PROGRAM FILTER 3 

DIMENSION Y{300)»G1(30G)»G2(300)»G3(300) , G4 ( 33 C ) , YERROR ( 300 ) , 

2X1 (300) ,X2( 300) ,X 3(300 ) ,XN( 300) , TIME (300) , I TITLE ( 1 2 ) , P 1 1 ( 300) , 

3 PI 2 ( 300 ), PI 3(300) , PI 4(300 ), P21 ( 300) ,P22( 300) , P23(300) ,P24 (300) , 
4P3 1 (300) , P32 ( 300 ) , P33( 3C0 ) , P34( 300) , P41 (300) ,P42(300) ,P43(300) , 
5P44( 300) ,XNT(300) , XVT (3C0) , X2T( 300) , X 1 E ( 300 ) , X2 E ( 300 ) , X 1 E S ( 3CC ) , 
6 AVX 1 E ( 300 ) , VAX1 ES( 300) , VAX! E ( 3C0) , AVXNTt 300 , VAXNT (300), 

7XNTS (300 ), VAXNT S( 300 ), SNOOT H( 300) ,YERRS( 300) ,A VYERR ( 300 ) , 
8VAYERRS ( 300 ) , VAYERR( 300 ) ,ABSXNT ( 300 ) , ABSX1 E (300 ) 



C 

C 



READS COE 
VARIANCE. 



WHICH IS THE SQUARE ROOT OF THE 



THE NEXT CARD 
MEASUREMENT NOISE 
READ 709, COE 
709 FORMAT (F10.0) 

READ 2, P 1 2 { 3 ) , PI 3 ( 3 ) , P 1 4 ( 3 ) , 

2 P2 1 ( 3 ) , P2 3(3) ,P24(3), 

3 P31 (3) ,P32( 3) , P 34 ( 3 ) , 

4 P41 (3) ,P42( 3) ,P43(3) 

2 FORMAT ( 1 2F5 . 0 ) 

READ 3,Q1 1 ,Q22,Q33,Q44 

3 FORMAT (4F20.0) 

THE NEXT CARD READS THE SAMPLING PERIOD. 

READ 5,T 

5 FORMAT (F10.0) 

THE NEXT CARD READS THE INDEXES WHICH LIMIT THE VARIOUS 
PIECES OF TRACK. 

READ 300, IMAX1 » IMIN2, I MAX 2, IM IN 3 , IM AX3 , I M I N4 , IMAX4 
300 .FORMAT (7110) r . 

THE NEXT CARD READS IN SUCCESSION THE INITIAL DISTANCE 
AND SPEED AND THE ACCEL L ERA T I ONS IN THE FOUR PIECES OF TRACK. 
READ 303,X1T0,X2T0,X3T0 1,X3T02, X3T03,X3T04 
303 FORMAT ( 6F 1 0 . 0 ) 

THE NEXT CARD READS THE RADIUS OF ROTATION AND THE RATE 
OF TURN FOR THE CIRCULAR PIECE GF TRACK. 

READ 705, R, ROT 
705 FORMAT (2F10.0) 

PRINT 999 



A“1 



Cont . 



/' 



999 FORMAT ( 1H1 ) 

PRINT 708 

708 FORMAT! /////10H COE // ) 

PRINT 707, COE 
707 FORMAT ( F 1 0 .2 ) 

PRINT 800 

800 FORMAT (/////7H T//) 

PRINT 36, T 
36 FORMAT (F10.2) 

PRINT 301 

301 FORMAT (/////71H IMAXl IMIN2 IMAX2 IMIN3 

13 IMIN4 IMAX4// ) 

PRINT 302, I MAX 1 , IMIN2, IMAX2, IMIN3»I MAX3, IMIN4, IMAX4 

302 FORMAT (7110) 

PRINT 304 

304 FORMAT (///// 72H X1T0 X2T0 X3TC1 X3TC2 

2 X3T03 X3T04 //) 

PRINT 305,XlT0,X2T0,X3TCl,X3T02,X3T03tX3T04 



305 


FORMAT 


(6F12.8) 








PRINT 


306 






306 


FORMAT 


( /////28H R 


ROT //) 






PRINT 


307, R, ROT 






307 


FORMAT 


(2F14.8) 








PRINT 


706 






706 


FORMAT 


( /////56H Q 1 1 


Q22 


Q33 


2C44 


//) 






PRINT 


31 ,Q1 1 ,Q22,Q33,Q44 






31 


FORMAT 


(4F14.8) 








DO 88 


K=1 ,250 








TIME ( 


K)=K 






88 


TIMEIK 


)=T*TIME(K) 







Cont . 



A”2 



I MAX 



ooon 



THE NEXT GROUP OF CARDS EVALUATES THE OUTPUT OF THE 
OBSERVABLE PROCESS. X1T AND X2T ARE THE REAL VALUES OF 
DISTANCE AND SPEED, XNT IS THE MEASUREMENT NOISE AND Y IS 
THE OUTPUT OF THE OBSERVABLE PROCESS. 

NUN IF=I22 07 03125 
DO 105 1=1 , I M AX 1 
CALL RNDEV (NUNIF.CEV) 

XNT ( I ) =COE*DEV 

Xl T ( I ) =X1T 0 +R*S IN F C ROT* ( T IME(I)-T I ME ( 1 ))) 

X2T ( I ) =X2T0 * COS F { ROT* ( T I ME ( I) — T IM E ( 1 )J) 

105 Y ( I ) =X 1 T { I ) +XNT ( I ) 

DO 102 1= IMIN2 , IMAX2 
CALL RNDEV { NUN IF , DEV ) 

XNT ( I ) =COE*DEV 

XI T( I ) = X 1 T ( I MAX 1 ) + { T I ME { I )-TIME( I MAXI ) ) *X2T( I MAXI ) + 

2 0.5*(TIME(I )-TIME( IMAX1) )**2*X3T02 

X2T( I ) = X2T ( I MAX 1 ) + ( T I ME ( I )-T IME ( I MAX 1 ) )*X3T02 

102 Y ( I ) = X 1 T { I ) +XNT ( I ) 

DO 103 1= IMIN3, IMAX3 
CALL RNDEVf NUNI F,DEV ) 

XNT { I )=COE*DEV 

XI T( I )=XlTt IMAX2)+(TIMECI)-TIME( IMAX2) ) *X2T( IMAX21+ 

2 0.5*{TIME( Il-TIMEI IMAX2 ) )**2*X3T03 

X2T( I ) = X2T ( IM AX'2 ) + ( T IME ( I I — T IME ( IMAX2) )*X3T03 

103 Y ( I ) = X 1 T ( I ) +XNT ( I ) 

DO 104 I = I M I N4 , IMAX4 
CALL RNDEV (NUNIF,OEV) 

XNT ( I )=COE*DEV 

XI T { I ) =X 1 T ( I MAX3I + (TIME (I l-TIMEC IMAX3) ) *,X2T ( IMAX3 ) + 

2 0.5*(TIME( I l-TIMEC IMAX3) )**2*X3T04 

X2T ( I ) =X2T( IMAX3) + (T IME ( I )- T I ME ( IMAX3 ) ) *X3T04 

104 Y { I ) =X 1 T ( D+XNTII) 

XI ( 3 ) =Y ( 3 ) 

X2 (3 ) = ( Y ( 3 ) — Y { 2 ) )/T 

X3(3)=( ( (Y(3) — YC2) )/T)-((Y( 2 ) — Y ( 1)) / T ) )/T 
XN ( 3 ) =0 . 0 
PI 1 ( 3 ) =Q44 

P22( 31=2.0/ {T**2)*Q44 
P33(3)=6.0/(T**4)*Q44 
P44 ( 3 ) =Q44 
X1ES(2)=0.0 
VAX 1 ES ( 2 ) =0 . 0 
XNrS(2)=0.0 
VAXNTS ( 2 ) =0 . 0 
YE RRS ( 2 ) =0 . 0 
VAYERRS ( 2)=0.0 
600 CONTINUE 
1 = 3 
J=I-1 
GO TO 10 
21 J= I - 1 



Cont . 



THE NEXT GROUP OF CARDS EVALUATES THE COVARIANCE MATRIX 
AJl^l.O-Gl ( J) ) 

Al3=ol6 

A1L=-G! ( J) 

A21=-G2(J) 

A22= 1 .0 
A 2 3 = 0 . 0 
A24=-G2(U) 

A3 1 =-G3 ( J ) 

A32=0.0 
A33= 1 .C 
A34=-G3 { J 1 
A4 1 =-G4 ( J ) 

A 4 2 = 0 . 0 
A43= 0 . 0 
A44= 1 . 0-G4 ( J ) 

D 1 1 = P 1 ] ( J)+P12{ J)*T+P13{ J )*0.5*T**2 
B12=P12(J)+P13{ J)*T 
B1 3=P 1 3 ( J ) 

614=0.0 

B2 1 = P21 (J)+P22(J)*T+P23(J )*0.5*T**2 

B22=P22(J)+P23(J)*T 

B23=P25( J) 

B24=0.0 

B31=P31 ( J)+P32( J )*T+P33( J )*0.5*T**2 
B32=P32{ J )+P33{ J)*T 
B33=P33 ( J ) 

834=0.0 

B4 I =P4 1 ( J )+P42{ J) *T + P43( J )*0.5*T**2 
B42=P42( J)+P43( J)*T 
B43=P43 ( J ) 
e44=0.0 

Cl 1 = A1 1 *B1 1 + A 1 4«-B*4 1 
C 1 2= A 1 1 *B 1 2+ A 1 4*842 
C13=A1 1*B13+A14*E43 
C 1 4=0 . 0 

C21=A21*B1 1+A22*B21+A24*B41 
C22=A21*B12+A22*B22+A24*B42 
C23=A21*B13+A22*B23+A24*B43 
C24=0 . 0 

C31=A31 *B1 1+A33*B31+A34«B41 
C32=A31*B12+A33*B32+A34*B42 
C33=A31*B13+A33*B33+A34*843 
C34=0.C 

C41=A41 *B1 1+A44*B41 
C42=A4 1*B12+A44*B42 
C43=A41 *B1 3+A44*e43 
C44=0.G 

PI 1 ( I )=C1 1+T*C21+0.5*T**2*C31+Q1 1 
PI 2 ( I )=C12+T*C22+0.5*T**2*C32 
PI 3 ( I )=C13+T*C23+0.5*T**2*C33 

P 1 4 ( I ) =0 . 0 r. 

P21 ( I ) = C2 1 +T*C3 1 

P22 { I )=C22+T*C32+Q22 

P231 I ) =C23+T *C33 

P24 ( I )=0.0 

P3 1 ( I ) =C3 1 

P32 ( I ) = C32 

P33 ( I ) = C33+Q33 

P34 ( I ) = 0 .0 

P41 ( I )=0.0 

P42 ( I ) =0 . 0 

P43 { I )=0.0 

P44 ( I ) = Q44 Cont. 



A =4 



on 



THE NEXT GROUP OF CARDS EVALUATES THE GAINS AND THE 
FILTERtD STATE VECTOR. 

10 S I GMAS = P 1 1 ( I ) + P 1 4 ( I ) + P4 1 ( I ) +P44 ( I ) 

G1 ( I ) = ( P 1 1 ( I ) + P14( I) )/S IGMAS 
G2 ( I )=( P21 ( I ) +P24 ( I ) ) /SIGMAS 
G3(I ) = ( P3 1 ( I ) + P34 ( I) )/S IGMAS 
G4( I )=(P41 ( I ) + P44( I) )/S IGMAS 
IF (1-3) 200,200,201 

200 YERRGR(I)=Y( I ) - ( X 1 { I )+XN( I) ) 

XI ( I ) = X1 ( I)+YERRCR(I ) *G1 ( I) 

X2 ( I ) = X2 ( I ) +YERROR ( I ) *G2 ( I ) 

X3(I ) = X3 ( I)+YERROR( I )*G3( I) 

XN( I )=XN( I ) +YERROR ( I )*G4( I ) 

GO TO 202 

201 YERROR( I ) =Y ( I ) - ( X 1 ( J ) +T *X2( J ) + 0 .5*T**2*.X3( J > ) 

XI ( I )=X1 ( J)+T*X2( J )+0.5*T**2«X3(J ) + YERROR( I ) *G 1 ( I ) 

X2( I ) = X2 ( J ) +T *X3 ( J ) + YERROR( I ) *G2 ( I ) 

X3( I )=X3( J )+YERROR( I ) »G2( I ) 

XN ( I ) =XN( J J + YERROR ( I )«G4( I) 

202 XI E( I )=X1 ( I ) -X 1 T t I ) 

X2E ( I ) =X2 ( I ) -X2T { I ) 

TOTI = I -2 

XI ES ( I ) = X 1 ES ( J) +X1 E( I ) 

AVX 1 E ( I ) =X 1 ES ( I J/TCTI 
VAX1ES(I)=VAX1ES(J ) + X 1 E ( I ) * *2 
VAX 1 E ( I ) = VAX 1 ES ( I ) /TOTI 
XNTS ( I )=XNTS ( J) +XNT( I ) 

AVXNT ( I )=XNTS( I )/TGTI 
VAXNTS ( I ) =VAXNTS ( J ) + XNT ( I ) * *2 
VAXNTf I ) = VAXNT S ( I ) /TOTI 
YERRS ( I ) =YERRS( J) +YERROR( I ) 

AVYERRf I )=YERRS ( I ) /TOTI 

VAYERRS( I )=VAYERRS ( J )+YERROR( I ) **2 

VA YE RR ( I ) =V AYERRS ( I ) /TOTI 

SMOOTH ( I ) = ( ( VAXNT ( I ) -VAX 1 E ( I ) ) /VAXNT ( I ) )*1 00.0 
ABSXNT ( I )=ABSF(XNT ( I ) )+6. 0 
ABSX1 E( I ) =-ABS F ( X 1 E ( I ) ) +6.0 
1 = 1 + 1 

IF ( I — I MAX4- 1 ) 21,22,2 2 
22 CONTINUE 
PRINT 999 



A~§ 
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THE NEXT GROUP OF CARDS PRINTS OUT THE RESULTS, 

PRINT 756 

756 FORMAT (////// 61H 

, 2 Q33/// ) 

\ PRINT 757, Q33 

757 FORMAT ( F68 .15) 

PRINT 751 

751 FORMAT (///119H TIME XNT XI T XI XI E 

2 X2T X2 X2E X3 AVXNT AVX1E SMCOT 

3H // ) 

PRINT 750, ( ( TIME( I ) , XNT ( I ) , X 1 T ( I ) ,X 1 ( I ) , X 1 E ( I ) ,X2T< I ) , X2 ( I ) , 

2 X2E ( I ) , X3 ( I), AVXNT (I ), AVX1 E( I ) , SMOOTH! I ) ) , 1 = 4, IMAX4) 

750 FORMAT ( F9. 1 , 1 1 E10.3 ) 

PRINT 999 
PRINT 752 

752 FORMAT (//////69H TIME P 1 1 P22 P33 G1 

2 G2 G3 //) 

PRINT 753, ((TIME(I),P11< I),P22( I ) ,P33 ( II ,G 1 (II , G2 ( I ) , G3 (I) ) , 

2 1=4, IM AX4 , 2 ) 

753 FORMAT (F9. 1 ,6E10.3) 

DO 11 L= 1 » 1 2 

11 IT IT LE ( L ) =8H 

I TITLE I 1 )=8HCARL0S P 
I T I T LE ( 2)=8H. SIMOES 
ITITLE(4)=8H MS T 
I T I TL E ( 5 J =8HHES I S 
I T I T LE( 7 ) =8H OP 

ITITLE(8)=8HTIMAL FI 
ITITLE(9)=8HLTERING 
ITITLEC 10) =8H0F RADAR 
ITITLE ( 1 1 ) =8H DATA 
LA6EL-4H P11 

CALL DRAW ( 1 49, T I ME , P 1 1 , 1 , 0 , L ABEL ,1 T I TLE, 0, 0, 0, 0, 0, 0 , 5 , 6, 1 , LAST ) 

LABEL=4H P22 

CALL DRAW ( 1 49 , TIME, P22 , 2 , 0 ,L ABEL ) 

LABE L=4H P33 

CALL DRAW ( 1 49, T IME , P33 , 3, 0 , LAB EL 1 
LABEL=4HVXNT 

CALL DRAW (249, TIME, VAXNT , 1 ,0, LABEL, I T I TLE , 0 , 0 , 0 , 0, 0 , 0 , 5 , 6 , 1 .LAST) 
LABEL=4HVX1 E 

CALL DRAW ( 2 49 , T IME , V AX IE , 2 , 0, LAB EL ) 

LABEL=4H 

CALL DRAW ( 249 , T I ME, AVXNT , 2 , 0 , L ABEL ) 

LABEL = 4H 

CALL DRAW ( 249, T I ME , AVX 1 E ,3 , 0 , L ABEL ) 

CALL DRAW(249,TIME , ABSXNT , 1 ,0, LABEL, I T I TL E , 0 ,0 , 0 , 0, 0, 0. 5 ,6 , 1 , LAST ) 
LABEL=4HX 1 E 

CALL DRAW (249, T IME, ABSX IE, 3,0, LABEL ) 

LABEL=4 HSMOO 

CALL DR AW (249, TIME, SMOOTH, 0,0, LABEL, ITI TLE , 0 ,0 , 0, 0 , 0, 0, 5 ,6 , 1 , LAST ) 
601 CONTINUE 
END 



INPUT DATA 



COE 

2.00 

T 

2.00 

I MAXI 
60 

XI TO 

7.00000000 

R 

6.40000000 
Q1 1 



IMIN2 

61 



I MAX 2 
148 



IMIN3 

149 



IMAX3 

162 



IMIN4 

163 



IMAX4 

250 



X2T0 

.27800000 



X3T01 


X3T02 


X3T03 


X3T04 


OOCOOOOO 


.00000000 


-.00983000 


.00000000 



RCT 

.05237000 

Q22 Q33 C44 



00000000 



00C00000 



0000 09 00 



4.0CCC0C00 



Case d) 
OUTPUT DATA 



TIME 


XNT 


X IT 


XI 


X 1 E 


8.0 


.263E+01 


. 898E + 0 1 


. 106E+02 


. 166E+01 


10.0 


-.11 5 E-0 1 


. 960E+0 1 


.971 E + 0 I 


. 103E + 00 


12.0 


-. 158E+00 


. 102E+02 


.976E+01 


-. 445E+00 


14.0 


.224E+01 


. 108E+02 


. 123E+02 


. 1 5 1 E+0 1 


16.0 


- . 2^7 E+0 1 


. 1 13E+02 


• 98 9 E+0 1 


-. 140E+01 


1 8.0 


-. 190E+01 


. 1 18E+02 


.958E+0 1 


-.213E+01 


20.0 


-. 1 1 9E+01 


. 122E+02 


. 102G + 02 


-. 195E+01 


22.0 


-.2C9E+Q1 


. 125E+02 


. 102E + 02 


-. 230E+0 1 


24.0 


-. 136E+01 


. 128E + 02 


. 109E+02 


-. 198E+01 


26.0 


-. 150E+01 


. 131E + 02 


. 113E+02 


-. 1 8 1 E+0 1 


28.0 


.84 2E+00 


. 133E + 02 


. 129E + 02 


-.4C6E+00 


30.0 


-.1 blE+01 


. 134C + 02 


. 127E+02 


— • 687E +00 


32.0 


. 394E+0 1 


. 134E + 02 


. 151E+02 


. 174E+01 


34.0 


- . 335E+00 


. 134E + 02 


. 148E+02 


. 141E+01 


36.0 


. 425E+0 1 


. 133E+02 


. 164E + 02 


.3152+01 


38.0 


-. 198E+0 1 


. 151E+02 


. 148E+02 


. 175E+01 


40.0 


. 8 1 6 E+00 


. 128E + 02 


. I 47 E+0 2 


. 1 81E+01 


42.0 


.418E+01 


. 125E+02 


. 157E+02 


. 312E+01 


44.0 


-. 152E+01 


. 122E+02 


. 14 1E+02 


. 1 95E+01 


46.0 


- . 1 90E+0 1 


. 1 1 8E+02 


• 1 26E+0 2 


. 884E+00 


48.0 


-. 106E+01 


. 1 13E+02 


. 1 16E+02 


. 332E+00 


50.0 


- . 3 1 6E+00 


. 108E+02 


.-109E + 02 


. 1 37E+00 


52.0 


.21 5E+01 


. 102E+02 


. 1 1 0E + 02 


. 783E+00 


54.0 


-.76 1E+00 


. 960 E + 0 1 


.993E+0 1 


. 327E+00 


56.0 


- . 926E+00 


. 897E+0 1 


.886E+01 


-. 1 15E + 00 


58.0 


.21 5E+01 


.833E+01 


. 880E+0 1 


. 478E+00 


60.0 


-.951 E+00 


. 767 E+0 I 


. 764E+0 1 


-. 29 1 8-01 


62.0 


.850E+00 


. 700 E + 0 1 


.709E+01 


. 924E-0 1 


64.0 


.38 5 E-0 1 


. 633E+0 1 


. 625E+0 1 


-. 765E-0 1 


66.0 


. 482E+00 


.567 E+0 1 


. 557E+0 1 


-. 977E-0 1 


68.0 


- • 1 3 1 E+0 1 


. 502E+0 1 


. 434E+0 1 


-. 679E+00 


70.0 


-.222E+01 


. 439 E+0 1 


. 294E+0 1 


-. 1 45 E + 0 1 


72.0 


. 4 1 3E+00 


. 380 E+0 1 


. 25 1 E +0 1 


-. 1 29E+01 


74.0 


-. 1 91 E+00 


. 323E+0 1 


. 191E+0 1 


-. 133E+01 


76.0 


-.1 84E-01 


. 27 1 E+0 1 


.14 1 E+0 1 


— . 1 3 1 E+0 1 


78.0 


-.624E+00 


. 224E + 0 1 


. 778E+00 


-. 146E + 01 


80.0 


-, 145E+01 


. 182E+01 


-.171 E-0 1 


-. 184E+01 


82.0 


.700E+00 


. 145 E+0 1 


-.214E-01 


-. 1482+01 


84.0 


-. 145E+00 


. 1 15E + 01 


-. 268E+00 


-. 1 42E + 01 


86.0 


-.171 E+0 1 


. 91 2E+00 


-.922E+00 


-. 183E+01 


88.0 


. 123E+01 


. 739E+00 


-. 52 7E + 00 


-.1272+01 


90.0 


. 262E+0 1 


. 634 E+00 


.286E+0C 


-. 348E+00 


92.0 


.2 1 5 E+00 


. 600 E + 00 


. 289E+00 


-. 3 11 E+00 


94.0 


. 168 E+0 1 


. 636E+00 


. 807E+00 


. 1 7 1 E+00 


96.0 


. 288 E+0 1 


. 74 1 E+00 


. 168E+0 1 


. 941 E+00 


98.0 


. 396E+0 1 


.915 E+00 


.282E+01 


. 1 91E+01 



A -8 



X2T 



.264E+00 
. 254E + 00 
.24 1E + 00 

* 225E +00 
•2C7E+00 
. 1 86E + 00 
. 163E + 00 
. 1 39E+00 
. 1 13E + 00 
. 858E-01 
.577E-01 
. 290E-0 1 
.844E-04 
.29 IE-01 
. 579E-0 1 
. 860E-0 1 
.1 13E+00 
. 139E+00 
. 163E+00 
. 186E + 00 
.2C7E+00 
. 225E+00 
.241E+00 
.2 54E + 00 
.264E+00 
.272E+00 
. 276E+00 
■.278E+00 
.276E+00 
. 272E + 00 
. 264E+00 
•.254E + 00 
.24 1E + 00 
• . 225E + 00 
■. 206E+00 
■. 1 8 6E + 00 
• . 1 63E+00 

139E+00 
•. 1 13E + 00 
• . 857E-0 1 

• .576E-01 
■.288E-01 

. 253E-03 
.293E-01 
.581 E-0 1 
.862E-01 



X2 



• 332E+00 
. 703E+CC 
.44 1 fc+CC 
. 662E+C0 
. 443E+C0 
.421E+C0 
. 155E+00 
. 129E+00 
. 146E-C1 
.761E-01 
. 3 1 4E+00 
. 218E+00 
.523E+00 
.379E+C0 
.5 15E+00 
. 239E+00 
. 173E+C0 
.246E+C0 
.408E-01 
. 128E+00 
•226E+CC 
.280E+C0 
. 250E+00 
. 327E+CC 
.395E+C0 
.360E+CC 
. 428E+CC 
. 430E+C0 
.455E+CC 
.461 E+CC 
. 5 1 4E+00 
. 575E+C0 
. 539E+C0 
. 520E+C0 
.491E+00 
. 475E+C0 
. 474E + C0 
.4 01 E + CO 
. 355E+00 
.351E+00 
.25 1E + C0 
. 124E + C0 
. 836E-0 1 
. 207E-02 
. 1 02E+CC 
. 2 1 7E+00 



X2E X3 AVXNT AVX 1 E SMOOTH 



.673E-01 
.957E+00 
.682E+00 
.457E+00 
.649E+00 
.607E+00 
. 3 1 9E+00 
.268E+00 
.984E-01 
. 975E-C2 
. 2 56 E + CO 
.189E+00 
.523E+00 
. 408E+00 
.573 E + CO 
.325E+00 
. 286E+00 
. 385E+00 
. 204 E+CO 
.585 E-G 1 
. 194E-01 
.552E-01 
.957E-C2 
.732E-01 
. 130E+00 
. 882E-0 1 
. 151E+00 
. 152E+C0 
. 179E+00 
. 190E+00 
. 250E + 00 
.321 E+00 
. 299 E+00 
.295E+00 
.285E+00 
.289 E+00 
.311 E+00 
.2 62 E+CO 
. 242 E+00 
.265 E + 00 
. 194E+00 
.953E-01 
•.839E-01 
.314E-01 
. 439E-01 
.1 31 E+00 



-. 1 96 E + CO 
-. 360E+C0 
-. 159E+00 
. 797E-0 1 
-.981 E-01 
-.717E-01 
-.287E-01 
-. 209E-C 1 
-.533E-C2 
. 165E-03 
. 167E-01 
.840E-02 
. 258 E-01 
. 147E-01 
. 204E-0 1 
.431E-02 
.653E-03 
•402E-02 
-. 572E-02 
-. 127E-01 
-. 160E-01 
-. 169E-01 
-. 142E-01 
-. 163E-01 
-. 178E-01 
-. 147E-C1 
-. 164E-C1 
-. 150E-01 
-. 148E-01 
-. 138E-C1 
-. 149E-01 
-. 163E-C1 
-. 132E-01 
-. 1 12E-C1 
-.883E-C2 
-.731 E-02 
-.661E-02 
-. 267 E-02 
-. 368E-03 
-. 138E-03 
.437 £-02 
.973 E-02 
. 107E-01 
. 134E-01 
. 169E-01 
. 206E-0 1 



. 170E+01 
. 113E+01 
.805E+00 
. 109E+01 
. 499 E+00 
. 157E+00 
-. 1 14E-01 
-.243E+00 
-. 355E+00 
-. 459E+00 
- . 350E+00 
- • 440 E + 00 
-. 127E+00 
-.141 E+00 
. 134E+00 
.9 32 E-02 
. 542E-01 
.271 E+00 
. 182E + 00 
, .826E-01 
• 305E-0 1 
. 154E-01 
. 104E+00 
.. 698 E-0 1 
.315E-01 
. 1 1 0 E+00 
.721 E-0 1 
.989 E-0 1 
. 969E-0 1 
. 109E+00 
.650E-01 
-.418 E-02 
. 809E-02 
.241 E-02 
. 183 E-02 
-. 151 E-01 
528E-01 
-.335E-01 
- .363E-01 
-.771 E-0 1 
- . 459E-0 1 
. 162 E-01 
. 207 E-01 
. 577E-0 1 
. 1 1 9E+00 
.201 E+00 



. 121E+C1 

.839E+0Q 
.51 8E+C0 
. 71 6E+00 
. 364E+00 
. 607E-03 
-. 243E+00 
-.472E+C0 
- .6 23E + 00 
-.730E+00 
-.703E+00 
- . 702E+00 

- .527E + 00 
-. 398E+00 
-. 176E+C0 
-.630E-01 

.413E-01 
.203E+00 
.291E+00 
.3 19E + 00 
.320E+00 
.312E+00 
.33 1E+00 
.33 1E+00 
.3 14E+00 
. 320E+C0 
.308E+C0 
.300E+00 
. 288E+00 
. 275E+C0 
•245E+00 
. 194E + 00 
. 1 50E + 00 
. 108E+00 
. 688E-0 1 
. 275E-0 1 
-.216E-01 
-. 589E-01 
-.929E-01 

- . 1 35E+00 
-. 162E+00 
-. 1 67E+00 
-. 170E + 00 
-. 162E+00 
-. 138E+00 
-.948E-01 



.558E+02 
.557E+02 
.532E+02 
. 538E+02 
.584E+C2 
.438E+02 
.31 1E+02 
• 229E+C2 
. 146E+02 
. 1C4E + 02 
. 1 18E+02 
. 1 62E + C2 
.359E+C2 
.321E+C2 
. 355E+02 
. 348E + 02 
.3C9E+02 
. 3 35E+02 
.31 1 E +02 
. 32 8E + C2 
.335E+02 
. 335E+02 
. 359E+02 
. 362E+C2 
. 367E+C2 
. 392E+02 
• 397E+02 
.40 1E + C2 
.4C1E+02 - 

.402E+02 
. 407E+02 
.41 4E + 02 
.40 1E+02 
. 386E+02 
.371E+02 
.355E+02 
. 338E + 0 2 
.323E+02 
. 306E + C2 
•295E+02 
.29 1E + 02 
. 327E+02 
.327E+C2 
.34 1E + 02 
.373E+C2 
.412E+02 



TIME 



XNT 



XI 



XI E 



X2T 



100. 


.0 


.274E+01 


102. 


.0 


.521 E+01 


1 04 . 


.0 • 


-.4C3E+01 


1 06. 


,0 


•253E+00 


108. 


.0 


•2C5E+00 


1 10. 


,0 


.17 1 E + C 1 


112. 


.0 ■ 


-.301E+01 


1 14. 


,0 • 


1 92E+0 1 


116. 


,0 


-.270E+01 


1 18. 


,0 • 


-. 150E+01 


120, 


.0 


- . 343E+00 


122. 


.0 


-.228E+G1 


124. 


.0 


1 54E+0 1 


126. 


.0 


.856E-01 


128. 


.0 • 


-.772E+00 


130. 


.0 • 


-. 1 63E+01 


1 32. 


.0 


. 982 E+00 


134. 


.0 


- . 398 E+G 1 


136. 


.0 • 


-.632E+00 


138. 


.0 


1 52 E+00 


140. 


.0 


-. 154E+01 


142. 


.0 • 


101 E+01 


144. 


.0 


168E+01 


146. 


.0 


.434E+01 


148. 


.0 


-.188E+01 


150, 


.0 


. 155 E+01 


152, 


.0 


. 103E+01 


154. 


.0 


.2S6E+01 


156. 


.0 


17 1E+01 


158. 


.0 


.349 E+00 


160, 


.0 ■ 


-. 166E+01 


162. 


.0 


-.2C3E+01 


164, 


.0 


• 5 32 E+00 


1 66. 


.0 


. 1 34E+0 1 


168. 


.0 


-.499E-01 


170. 


.0 


-. 1 44E+00 


172. 


.0 


- . 247E+0 1 


174. 


.0 


- . 449 E+0 1 


176. 


.0 


.143E+01 


178. 


.0 


.257E+01 


180. 


.0 


.250E+01 


1 82, 


.0 


472E+01 


184, 


.0 


. 1 62E+01 


186. 


.0 


. 1 09E+01 


188, 


.0 


- . 983E+00 


190. 


.0 ' 


-. 1 62E+00 


192. 


.0 


. 847E+00 


194, 


.0 


-.378E+01 


196, 


.0 


-. 163E+00 


198, 


.0 


• 209E+0 1 



X 1 T 

. 116E+01 
. 146E+01 
. 183E+01 
. 225E+0 1 
. 272E+01 
. 324E+01 
.381 E+01 
.440E+01 
. 503E+01 
. 568E+01 
.634E+01 
.689E+01 
.744E+01 
. 800E+0 1 
.855E+01 
.910 E+01 
.966E+01 
. 102E+02 
. 108E+02 
. 113E+02 
. 1 19E+02 
.124E+02 
. 130E+02 
. 135E+02 
. 1 41 E+02 
. 146E+02 
. 152E+02 
. 157E+02 
. 1 63E+02 
. 168E+02 
. 174E+02 
. 1 80 E + 02 
. 185E+02 
. 191 E+02 
. 196E+02 
. 202E+02 
. 207E+02 
.21 3E+02 
.218 E+02 
. 224E+02 
. 229E+02 
. 235 E+02 
.240E+C2 
. 246E+02 
.25 1 E+02 
. 257 E+02 
.262E+C2 
. 268E+02 
. 274 E+02 
.279E+02 



.348E+0 1 
.48 6 E+G 1 
. 329E+0 1 
. 337E + 0 1 
. 356E+0 1 
. 430E+0 1 
. 362E + 0 1 
•354E+C1 
. 339E + 0 1 
. 38 1 E +0 1 
. 469E + 0 I 
. 497E+0 1 
. 556E + 0 1 
. 667E+0 1 
.745E+0I 
. 795E+0 1 
. 926E+0 1 
. 893E+0 1 
• 977E + 0 1 
. 107E+02 
. 1 1 2E+02 
. 1 1 8E+02 
. 122E+02 
. 144E+02 
. 145E+02 
. 156E+02 
. 165E+02 
. 177E+02 
. 175E+02 
. 180E+02 
. 179E+02 
. 177E+02 
. 185E+02 
. 194E+02 
. 198E+02 
. 203E+02 
. 200E+02 
. 193E+02 
.206E+02 
.22 1E + 02 
. 235E+02 
.225E+02 
.237E+02 
. 247E+02 
. 249E+02 
.255E+02 
.264E+02 
. 258E+02 
. 265E+02 
. 27 8E + 02 



.232E+01 
. 340E+01 
. 146E+01 
. 1 1 3E+01 
. 834E + 00 
. 1 C6E+01 
-. 187E+00 
-. 868E + 00 
-. 1 64E+01 
-. 187E+01 
-. 1 65E+01 
-. 192E+01 
-. 1 89E+01 
-. 1 33E+01 
-. 1 10E+01 
-. 1 1 6E+01 
4C0E+00 
-. 1 28E+01 
-. 992E+00 
-. 605E+00 
-.71 0E+00 
-. 648E+00 
-. 808E+00 
. 884E+00 
. 377E+00 
. 969E+00 
. 128E+01 
. 1 97E+01 
. 1 20E+01 
. 1 15E+01 
. 468E+00 
-.237E+00 
-.54 IE-01 
. 339E+00 
. 236E+00 
. 1 1 3E+00 
-. 7C4E+00 
-. 200E+01 
-. 125E+01 
-. 264E+C0 
.551 E+00 
-. 974E + 00 
-. 309E+00 
. 893E-0 1 
-. 208E+00 
-. 206E+00 
. 105E+00 
-. 1 05E+01 
-. 899E+00 
-. 776E-0 1 



71T3E + 00 
. 1 39E+00 
. 1 64E+00 
. 186E+00 
. 207E+00 
.225E+00 
.24 1 E + 00 
. 254E+00 
. 264E+00 
.272E+00 
. 277E+00 
.277E+00 
.277E+00 
.277E+00 
. 277E+ 00 
. 277E+00 
. 277E+ 00 
.277E+00 
.277E+00 
. 277E+00 
.277E+00 
.277E+00 
• 277E + 00 
.277E+00 
. 277E+00 
.277E+00 
. 277E + 00 
. 277E+00 
. 277E+ 00 
. 277E+00 
. 277E + 00 
. 277E+00 
.277E+00 
. 277E+00 
. 277 E+00 
. 277E+00 
.277E+00 
. 277E+00 
. 277E + 00 
. 277E+00 
.277E+00 
. 277E+00 
. 277E+00 
.277E+00 
.277E+00 
. 277E+00 
.277E+00 
. 277E + 00 
.277E+00 
. 277E+00 



A»10 



X2 

T275E+00 
. 390E + C0 
.22 1E + C0 
.216E+00 
.217E+00 
.2676+CC 
. 1 81E + 0C 
. 155E + 00 
. 125E+CC 
. 1 48E+C0 
.21 OE+CC 
.21 1E+CO 
. 238E+CC 
. 309E+00 
.341 E+CC 
. 343E+00 
.41 8E+CC 
. 334E + C0 
. 363E+00 
• 397E+00 
. 381 E+OC 
.381E+0G 
. 359E+00 
. 506E+00 
.439E+00 
.476E+CO 
.483E+00 
.522E+00 
.424E+00 
•399E+C0 
.31 9E+00 
. 244E+00 
. 258E+00 
. 290E+C0 
.274E+00 
• 259E+00 
. 1 82E+00 
. 708E-G 1 
. 156E + 00 
. 257E+00 
.333E+00 
. 1 89E + C0 
. 258E+00 
. 296E+C0 
. 268E+C0 
. 270E+CC 
. 299E+00 
. 193E+00 
• 21 5E + C0 
.296E+00 



X2E 



X3 



AVXNT 



AVX 1 E 



SMGOTH 



. 162E+00 
.250E+00 
.572E-01 
. 293E-01 
. 1 01 E-01 
.417E-01 
-. 601 E-01 
-.989E-01 
-. 1 39E+00 
-. 124E+00 
-. 664E-0 1 
-.659E-0 1 
-.388 E-01 
.320 E-01 
. 648E-0 1 
.670E-01 
. 141 E+00 
.577E-01 

• 864 E-0 1 
. 1 20E+00 
. 104 E+00 
. 1 04 E+00 

• 827 E-0 1 
. 230E+00 
. 163E+00 

• 200 E + 00 
. 206E+00 
.245 E+00 
. 147 E+00 
. 123 E+00 
. 425 E-01 

-.322 E-01 
-. 188 E-01 
. 131 E-01 
-. 246E-02 
-. 179 E-01 
-.942 E-0 1 
-. 206E + 00 
-. 120E+00 
-. 198 E-01 
.561 E-01 
-.873 E-0 1 
-. 186 E-01 
. 199 E-01 
-.841 E-02 
-. 698E-02 
.223E-01 
-.838 E-01 
-.619 E-0 1 
. 1 95E-01 



.213E-01 
. 246 E-01 
.147 £-0 1 
. 132E-01 
. 120E-01 
. 1 32E-0 1 
.8 10E-02 
.621 E-02 
. 429E-02 
.494 E-02 
. 729 E-02 
. 666E-02 
• 728E-C2 
. 983 E-02 
. 104E-01 
. 958E-C2 
. 121 E-Cl 
. 7 20 E-02 
. 784 E-02 
. 866E-02 
. 7 1 7E-C2 
.651 E-02 
. 494 E-02 
. 1 12E-C1 
.71 IE-02 
.8 15E-02 
.770 E-02 
. 877E-02 
.354 E-02 
.2 10E-02 
-. 172E-02 
-.495E-C2 
-. 390E-02 
-.2 10E-02 
-.261 E-02 
-. 307E-02 
-.625 E-02 
-. 1C7E-01 
-. 590E-02 
-. 8C8E-03 
.271 E-02 
-. 404 E-02 
-. 557E-03 
. 123E-C2 
-. 159E-03 
-.798E-04 
. 1 25 E-02 
-. 367 E-02 
-.234 E-02 
. 156E-02 



• 254E+00 
. 355E+00 
. 267E+00 
. 267 E+00 
. 266E+00 
. 293E+00 
.232E+00 
. 192E+00 
. 141 E+00 
. 1 12E+00 

• 1 04 E + 00 
.637 E-01 

• 370 E-0 1 
. 378E-0 1 
.247E-01 

-. 152E-02 
. 1 39 E-0 1 
-. 4 77E- 0 1 
-.56 5 E-01 

■ .580E-0 1 

797E-01 
-. 933E-01 
-. 1 16E + 00 
-.531 E-01 
-. 785E-0 1 
-. 562E-01 
-.416 E— 0 1 
-. 685E-02 
-. 293 E-01 
-. 244E-0 1 
- . 453E-0 1 
-.704E-01 
629 E-0 1 
-. 456E-01 
-.456E-01 
-.468E-01 
757E-0 1 
-. 128E+00 
-. 109E+00 
■ . 787E-0 1 
-.493E-0 1 
-. 1 02E+00 

■ .826 E-0 1 
• .698 E-0 1 
-. 797E-0 1 
■ . 806E-0 1 
-.707 E-0 1 
-. 1 10E+00 

110E+00 
-. 877E-0 1 



• .444E-0 1 
. 259E-0 1 
. 546E-0 1 
.756E-01 
. 902E-0 1 
. 108E+00 
. 103E+00 
♦ 853E-0 1 
. 545 E-0 1 
. 208E-0 1 
-. 795E-C2 
-.404E-01 
-.712E-01 
• . 9 1 9E-0 1 
-. 108E + 00 
-. 125E + 00 
•. 129E + 00 
•. 147E+00 
-. 1 60 E + 00 
-.1 66E+00 
•. 174E+00 
-. 181E + 00 
-. 190E+00 
-. 175E+00 
-. 167E+00 
-. 1 52E+00 
-. 1 32E + 00 
-. 104E+00 
• . 873E-0 1 
-.71 3E-0 1 

■ • 643E-0 1 
-. 665 E-0 1 
-.664E-01 
-.614E-01 
-. 577E-0 1 
-.557E-01 

■ .6 34E-0 1 
-.862E-01 
-. 997E-0 1 
-. 1 02 E + 00 
■ . 942E-0 1 
-.104E+C0 
-. 106E + 00 
-. 1 04E+00 
-. 105E + 00 
-. 106E + 00 
-. 1 04E+00 
-. 1 14E+00 
-. 122E+00 
-. 122E+00 



TT4 06E+C2 
.430E+02 
•463E+02 
. 458E+02 
. 454E + C2 
.457E+C2 
.479E+02 

• 484E + 02 
. 489E + C2 
.479E+C2 
. 467E + 02 
. 463E+02 
.454E+02 
.447E+02 
.443E+02 
. 444E+0 2 
.445E+02 
.472E4C2 
•470E+02 
. 468E+C2 
.471E+C2 
.472E+02 
•475E+02 
.506E+C2 
.5 1 2E+02 
.5 13E + 02 
. 509E + 02 
. 507E+C2 
. 507E + 0 2 
.5C3E+02 
.506E+02 
.51 3E+02 
.513E+02 
. 5 1 6E+02 
.51 5E+02 
. 5 1 5E+02 
•523E+02 
. 540E+02 
. 538E+02 
. 546E+02 
.553E+C2 
. 577E+02 

• 580E + 02 
.581E+02 
.582E+C2 
.582E+02 
.583E+02 
. 595E + 02 
.593E+C2 
•598E+02 



A=1 1 



TIME 



XNT 



XI 



X I E 



X2T 



200.0 -.801 E +00 

202.0 -.194E+00 

204.0 .181E+00 

206.0 -.269E+01 

208.0 . 1 T4E+0 1 

210.0 -. 166E+0I 

212.0 .197E+01 

214.0 .21 7E+01 

216.0 .11 9 E+0 1 

218.0 -.252E+00 

220.0 -.407E-01 

222.0 -.61 8E+00 

224.0 .4176+00 

226.0 .2 1 8 E+0 1 

228.0 . 254E+0 1 

230.0 -.357E-01 

232.0 .324E+00 

234.0 -.771E+00 

236.0 — . 293E+C 1 

238.0 .226E-C2 

240.0 -.137E+00 

242.0 -.274E+01 

244.0 . 1 57E+00 

246.0 .31 4E+01 

248.0 . 295E+0 1 

250.0 -.186 E+00 

252.0 . 427E+00 

254.0 . 1 54 E+0 1 

256.0 . 1 I 3 E+0 1 

258.0 .653E+00 

260.0 .1556+00 

262.0 . 120E+00 

264.0 .307 E+0 1 

266.0 -.61 3E+00 

268.0 - . 354 E+00 

270.0 -.250E+01 

272.0 . 160E+01 

274.0 . 1 46 E+0 1 

276.0 .71 2E+00 

278.0 -. 1 85E+0 1 

280.0 - .30 1 E+00 

282.0 . 254E+0 1 

284.0 -.373E+01 

286.0 - . 1 5 1 E+0 1 

288.0 .1426+00 

290.0 .7G0E+00 

292.0 .326E+01 

294.0 -.9446+00 

296.0 .26 2E + 00 

298.0 -. 192E+01 



X IT 

. 28 5E+02 
.290E+02 
. 296E + 02 
.301 E+02 
.3076+02 
.31 2E+02 
. 3 1 8E+02 
. 323 E+02 
.3296+02 
. 334 E+02 
. 340E+02 
.3456+02 
. 351 E+02 
. 356E+02 
.3626+02 
. 3686+02 
. 373E+02 
. 379E+02 
. 384E+02 
. 390 £+02 
. 395E+02 
.401 E+02 
. 406 E+0 2 
.4 1 2E+C2 
. 4 1 7 E+02 
. 423E+02 
.428E+02 
. 434 E+02 
.4396+02 
•445E+02 
.45 1E+02 
.4566+02 
. 462 E + 02 
.467E+02 
.473E+02 
. 478E+02 
.4846+02 
. 489E+02 
.4956+02 
. 500 E + 02 
. 506E+02 
.51 1 E+02 
.517 E + 02 
. 522 E + 02 
. 5286+02 
.5336+02 
. 539 E+02 
. 545E+02 
. 550E+02 
. 555 E+02 



.282E+C2 
. 288E+02 
.2956+02 
. 292E+02 
.3056+02 
. 306E+02 
.3196+02 
. 332E+02 
.3406+02 
. 343E+02 
. 347E+02 
.3496+02 
. 35 5E+02 
. 367E+02 
.3786+02 
.380E+02 
.384E+02 
. 38 4E+02 
.379E+02 
.384E+02 
.390E+02 
. 387E+02 
. 395E+02 
,.41 2E+02 
. 426E+02 
. 429E+02 
. 435E+02 
.4446+02 
. 45 1 E+02 
. 456E+02 
. 460E+02 
.4636+02 
. 47 6E+02 
. 477E+02 
.4796+02 
.47 5E+02 
. 48 5E+02 
. 494E+02 
. 500E+02 
. 498E+02 
.5026+02 
.51 5E+02 
.508E+02 
. 5 1 0E+02 
. 5 1 8 E +0 2 
.5276+02 
.544E+02 
. 546E+02 
. 55 2E+02 
. 552E+02 



-. 269E+00 
’-. 243E + 00 
-. 1 06 E+00 
-. 8696+00 
-. 1 49E+00 
-. 58 1 E + 00 
. 1656+00 
. 849E+00 
. 1 1 1 E + 0 1 
. 8696+00 
. 71 5E+00 
. 3986+00 
.4396+00 
. 1C0E + 01 
. 1 57E+01 
. 1256+01 
. 1066+01 
.5666+00 
-. 5 1 36+00 
-.5186+00 
-. 55 1 E+00 
-. 1 366+01 
-. 1 1 3E+01 
-. 9776-02 
. 878E+00 
. 662E+00 
. 661E+00 
.9956+00 
. 1 1 4 E + 0 1 
. 1 10E+01 
. 9106+00 
. 7296+00 
. 147E+01 
. 955E+00 
. 594E+00 
-.3656+00 
.8856-01 
. 4 33E + 00 
. 4976+00 
-. 2306+00 
-.3576+00 
. 415E+00 
-. 846E+00 
-. 121E+01 
-. 993E+00 
-.620E+00 
. 488E+00 
. 134E+00 
. 206E+00 
-. 387E+00 



. 277E+00 

.277E+00 

.277E+00 

. 277E+00 

.277E+00 

.2 77E + 00 

. 277E+00 

. 277E+00 

.277E+00 

. 2 77E + 00 

.2776+00 

.277E+00 

.277E+00 

.2 77 E+00 

•277E+00 

. 277E+00 

. 277E+00 

.277E+00 

.2776+00 

.277E+00 

. 277E+00 

.2776+00 

.2776+00 

.277E+00 

.2776+00 

.277E+00 

. 277E + 00 

. 277E + 00 

.277E+00 

.277E + 00 ' 

.277E+00 

.277E+00 

.277E+00 

.2776+00 

.2776+00 

.277E+00 

. 277E + 00 

.277E+00 

. 277E+00 

. 277E+00 

.277E+00 

.277E+00 

. 277 E+00 

. 277E+00 

.277E+00 

. 277E+00 

. 277E+00 

. 277E+00 

.2776+00 

. 257E+00 



A~12 



X2 



X3 



AVXNT 



AVX1E 



SMOOTH 



. 278E+00 
.281 E+CO 
. 294E + 00 
. 224E+00 
• 295E+00 
. 255E+00 
. 325E + CC 
. 383E+00 
. 397E+00 
. 362E+00 
•339E+CC 
. 302E + 00 
. 301 E+OO 
.348E+C0 
.390E+C0 

• 347E + C0 
. 320E+00 
.267E+C0 
. 1 67E+CC 
. 175E + C0 
. 1 80E+00 
. 115E + 00 
. 152E+C0 
. 267E+00 
. 350E+00 
.325E+00 
. 32 1 E+CO 
. 348E+00 
. 354E+00 
. 343E+00 
. 3 1 8E+00 
. 296E+00 
. 360E+00 
•303fc+00 
.266E+00 
. 1 78 E+OO 
.227E+00 
•263E+C0 
.269E+00 
. 204E+00 
. 199E + 00 
. 2 77E + C0 
. 164E + 00 
. 143E + 00 
. 1 77E+00 

• 223E + 00 
•332E+00 
. 298E + 00 
• 304E+00 
•248E+C0 



X2E 



.145E-02 
.460E-02 
. 1 74 E-0 1 
, 527E-01 
, 1 84 E-0 1 
.217 E-0 1 
. 487 E-0 1 
, 1 07E+00 
. 120E+00 
.859 E-0 1 
. 622 E-0 1 
. 2 55 E-0 1 
,245 E-0 1 
.710 E-0 1 
, 1 13E+00 
.701 E-0 1 
.437 E-0 1 
.908 E-02 
. 109E+00 
, 1 02 E+OO 
.968 E-0 1 
.161 E+OO 
. 125E + 00 
,994 E-02 
.734 E-0 1 
. 483 E-0 1 
, 446 E-0 1 
,711 E-0 1 
.775 E-0 1 
.662 E-0 1 
. 4 1 2 E-0 1 
. 194 E-0 1 
.832 E-0 1 
. 2 68 E-0 1 
. 1 05 E-0 1 
.983 E-0 1 
.493 E-0 1 
. 140E-01 
.719 E-02 
.726E-01 
.772 E-0 1 
.792E-03 
. 113 E + OO 
. 134E+00 
.992 E-01 
.531 E-01 
. 555E-0 1 
.211 E-01 
.279 E-01 
.872E-02 



cqq - 

» V I 4 - O 

.688E-C3 - 
. 121 E-02 - 
. 208 E-02 - 
. 1 33 E-02 - 
.606E-03 - 
. 264 E-02 - 
. 502E-02 - 
. 5 1 7E-02 - 
.316 E-02 - 
. 180E-02 ■ 
. 284E-04 - 
.686E-04 - 
.204E-G2 - 
.378E-C2 - 
. 148E-02 - 
.146E-C3 - 
.226 E-02 - 
. 660E-02 ■ 
. 566 E-02 - 
. 492 E-02 - 
. 739 E-02 - 
. 508 E-02 - 
.593E-03 - 
.431 E-02 ■ 
. 279E-C2 - 
. 237 E-02 - 
. 335 E-02 - 
. 334 E-02 
. 252 E-02 
. 1 17 £-02 
. 707E-04 
.295 E-02 
. 133E-03 
. 157E-02 
. 541 E-02 
.269E-02 
. 853E-03 
.467E-03 
.339 E-02 
.329 E-02 
.542E-03 
. 465E-02 
. 5 1 9E-02 ■ 
.315 E-02 
.771E-03 
. 422 E-02 
.228E-02 
. 238E-02 
.386E-03 



Q^np-n \ 

• W U V • 

. 960 E-0 1 
. 932E-0 1 
.11 9E+00 
. 101 E+00 
. 1 16E + 00 
. 957E-G 1 
.741 E-01 
.621 E-01 
. 639E-01 
■. 637E-0 1 
.6 88 E-0 1 
. 644E-0 1 
.442E-01 
.21 1 E-0 1 
.21 3 E-0 1 
. 1 82 E-0 1 
. 248E-0 1 
.498E-01 
.493E-01 
.501 E-01 
. 726E-0 1 
■ .707 E-0 1 
.442E-01 
•. 197E-01 
-.210E-01 
■. 174E-01 
•.492E-02 
.411 E-02 
• 922E-02 
.104 E-01 
.1 12E-01 
• 348E-01 
. 298E-0 1 
. 2 6 9 E- 0 1 
. 792 E-02 
. 198E-01 
. 304E-0 1 
.3 54 E-01 
.216 E-0 1 
. 193E-01 
. 375E-0 1 
. 106 E-01 
■.251 E-03 
. 748E-03 

• 564E-02 
.282 E-01 
.215E-01 
•232E-01 
. 994 E-02 



-s 123E+00 
-. 125E+00 
-. 124E + 00 
-. 132E+00 

- . 1 32E+00 
-. 136E+00 
-. 133E + 00 
-. 1 24E+00 
-. 1 12E + 00 
-. 103E+00 
-.956E-01 
-.91 IE-01 
-.862E-01 
-.764E-01 
-.617E-01 
-.50 IE-01 

- . 4 04 E-0 1 
-.351 E-0 1 
- • 392E-G 1 
-.433E-01 

- .476 E-0 1 
-. 586E-01 
-. 676E-01 
-. 67 1 E-0 1 
-.593E-01 
-.535E-01 
-.477E-01 
-. 394E-01 
- • 300 E-0 1 
-.211 E-0 1 
-. 1 38E-0 1 
-.805E-G2 

•332E-02 
. 106E-01 
. 1 50E-01 
.121 E-01 
. 127E-01 
. 158E-01 
. 1 94E-0 1 
. 176 E-01 
. 148 E-01 
. 177E-01 
.1 15E-01 
.286 E-02 
- .4 15E-02 
-.846E-02 
-. 50 1 E- 02 
- .405E-02 
-.261 E-02 
- . 523E-02 



598E+C? 
598E+02 
598E+C2 
603E+02 
606E+C2 
608E+02 
61 1E + C2 
6 1 4TE + 0 2 
612E+C2 
61 1E + 02 
61CE+C2 
610E+02 
609E+C2 
61 1E+C2 
61 1E + 02 
608E+02 
605E+C2 
605E+02 
612E+02 
61 1E + 02 
61 1E + 02 
6 13E+C2 
610E+02 
619E+02 
624E+C2 
623E+02 
622E+02 
622E+02 
62CE+02 
618E+C2 
616E+02 
615E+C2 
618E+02 
617E+C2 
616E+02 
621E+C2 
623E+02 
624E+02 
624E+02 
626E+C2 
626E+02 
631E+02 
639E+02 
638E+C2 
636E+C2 
636E+02 
642E+02 
643E+C2 
643E+C2 
645E+02 



A-13 



TIME 



300.0 

302.0 

304.0 

306.0 

308.0 

310.0 

312.0 

314.0 

316.0 

318.0 

320.0 

322.0 

324.0 

326.0 

328.0 

330.0 

332.0 

334.0 

336.0 

338.0 

340.0 

342.0 

344.0 

346.0 

348.0 

350.0 

352.0 

354.0 

356.0 

358.0 

360.0 

362.0 

364.0 

366.0 

368.0 

370.0 

372.0 

374.0 

376.0 

378.0 

380.0 

382.0 

384.0 

386.0 

388.0 

390.0 

392.0 

394.0 

396.0 

398.0 



XNT X IT 



XI X 1 E X2T 



. 1 45E+01 
.7 06E+00 
. 632E+00 
. 808E+00 
. 138E+01 
.1G1E+01 
. 141 E+00 
. 540E+00 
. 6580+00 
. 144E+01 
.21 2E+00 
. 985E+00 
. 937E-0 1 
.153E+01 
. 203E+0 1 
. 176E+01 
.349E+01 
.245E+01 
.688 E+00 
. 126E+01 
.248E+01 
.44 3 E+00 
. 144E+01 
.260E+00 
.7800+00 
•429E+0 1 
.6856+00 
.3300+01 
. 106E+01 
. 168E+01 
. 120E+01 
. 436E+00 
.338E+01 
. 175E+01 
. 175E+01 
.234E+01 
.1 05E+01 
. 104E+01 
.728E+00 
.223E+01 
.587E+00 
.779 E+00 
. 197E-01 
« 1 35E+01 
.226E+01 
.3 I7E+01 
. 974E+C0 
.142E+01 
. 275E+0 1 
.604E+00 



. 560 £+02 
. 565E+02 
. 569E+02 
.5730+02 
.576E+02 
. 579E+02 
.582E+02 
.584E+02 
. 586E+C2 
.5870+02 
. 588E+02 
.589E+02 
. 589E+02 
. 589E+02 
. 589E+02 
.5890+02 
. 589E+02 
. 589E+02 
. 589E+02 
. 589E+02 
. 589E+02 
. 589E+02 
. 589E+02 
.589E+02 
. 589E+02 
. 589E+02 
♦589E+02 
. 589E+02 
. 589E+02 
. 589E+02 
.589E+02 
. 589E+02 
. 589E+02 
. 589E+02 
.5890+02 
.590E+C2 
.59CE+02 
. 590E+02 
. 590E+02 
. 590E+02 
. 590E+02 
.590E+02 
.5906+02 
.5900+02 
. 590E+02 
. 590E+02 
.590E+02 
. 590E+02 
. 590E+02 
. 590E+02 



.55 3E+02 
. 558E + 02 
. 566E+02 
. 569E+02 
. 570E+0 2 
. 57 9E+02 
.582E+02 
. 588E+02 
. 592E + C2 
. 598E+02 
. 598E+02 
.601E+02 
. 599E+0 2 
.593E+02 
. 599E+02 
. 602E+02 
.610E+02 
.5980+02 
.5930+02 
.594E+02 
. 58 4E+02 
.5820+02 
. 586E+02 
. 586E+02 
.5870+02 
.6000+02 
.5990+02 
.586E+02 
.582E+02 
.5770+02 
. 582E+02 
.581E+02 
.593E+02 
.5970+02 
. 590E+02 
.5830+02 
. 587E+02 
. 584E+02 
.5876+02 
.581E+02 
.58 1E+02 
. 580E+02 
. 582E+02 
. 580E+02 
. 58 9E+02 
. 580E+02 
.5860+02 
.592E+02 
.60 1E+02 
.60 1E+02 



-.7100+00 
-.7250+00 
-. 303E+00 
-. 3720+00 
-. 5930+00 
-. 3980-01 
.7750-01 
. 3800+00 
.661 E+00 
. 1 120+01 
.9810+00 
. 1 210+01 
. 1 050+01 
. 436E+00 
. 971C+00 
. 1 30E+01 
.2060+01 
.8530+00 
. 3500+00 
. 510E+00 
— . 5 066 +00 
-.7296+00 
-. 3300+00 
-. 3446+00 
-. 1 85E+00 
. 1C3E+01 
.9450+00 
-. 3430+00 
-. 7310+00 
-. 1 23E+01 
-.7360+00 
-. 8040+00 
. 323E+00 
. 776E+00 
.8570-01 
-. 6696+00 
— . 2556+00 
-. 536E+00 
-. 222E+00 
-.8530+00 
-. 869E+00 
-. 927E+00 
-.7310+00 
-. 962E+00 
-.4520-01 
-.9280+00 
-. 400E+00 
. 182E+00 
. 1 07E+01 
. 1 14E+01 



.237E+00 
.21 8E+00 
.1 98E+00 
. 178E+00 
. 1 59E+00 
. 1 39E+00 
. 1 19E + 00 
.996E-01 
. 799E-0 1 
. 6C3E-0 1 
.4C6E-01 
. 2C9E-0 1 
. 1 27E-02 
. 1 27E-02 
. 127E-02 
. 127E-02 
. 127E-02 
. 127E-02 
. 1 27E-02 
. 127E-02 
. 127E-02 
. 1 27E-02 
. 1 27E-02 
.1270-02 
. 1 27E-02 
. 1 27E-02 
. 127E-02 
. 127E-02 
.1270-02 
.1270-02 
. 1 27E-02 
. 127E-02 
. 1 27E-02 
. 1 27E-02 
. 127E-02 
. 1 27E-02 
. 127E-02 
. 1 27E-02 
. 127E-02 
.1270-02 
. 1 27E-02 
• 1 27E-02 
.1270-02 
. 1 27E-02 
.1 27E-02 
. 127E-02 
. 1 27E-02 
. 1 27E-02 
. 1 27E-02 
. 127E-02 



A~!4 



X2 



. 2 1 8E+00 
.21 5E+C0 
.249E+00 
. 232E+C0 
. 199E+00 
• 236E+00 
. 22 7E+00 
. 23 1 E+ CO 
. 230E+00 
. 242E+00 
. 1 95E+00 
. 182E+00 
. 132E+00 
.450E-01 
•710E-01 
.770E-01 
. 123E+00 
. 1 38E-01 
. 723E-C1 
. 637E-0 1 
. 16 lE+OO 
. 175E+00 
. 129E+C0 
. 123E+C0 
. 101 E+00 
. 1 65E-01 
. 550E-C2 
.1 13E+00 
. 139E+00 
. 1 71 E+00 
.109E+CC 
. 102E+C0 
. 123E-C1 
.550E-01 
. 102E-C1 
. 756E-01 
. 287E-0 1 
. 489E-0 1 
. 132E-01 
. 666E-0 1 
. 589E-01 
. 557E-C 1 
. 295E-01 
* 444E— 0 1 
. 462E-0 1 
. 353E-C 1 
. 186E-01 
.721 E-01 
. 148E+00 
. 1 4 1E + CC 



X2E 



X3 



AVXNT 



AVX1E 



SMCOTH 



.1 92E-01 
.229 E-02 
.511 E-01 
. 535E-0 1 
.402E-01 
.969 E-01 
. 107E+00 
.152E + C0 
.151 E+00 
. 182E+00 
. 1 54 E+00 
, 1 61 E+00 
.131 £+00 
.437 E-01 
.698E-01 
. 757E-01 
. 122 E+00 
.150E-01 
. 7 36E-0 1 
.650 E-01 
. 162E+00 
. 176 E+OO 
. 1 30 E + 00 
. 124E+00 
. 1 02 E+00 
. 1 52 E-01 
. 422E-02 
. 11 5 E + OO 
. 140 E+00 
. 172E+00 
. 1 10 E+00 
. 1 04 E+00 
. 1 10 E-01 
.538E-01 
.115 E-01 
.768E-01 
.299 E-01 
.501 E-01 
.145E-01 
.678E-01 
.6 02 E-01 
.569E-01 
.307E-01 
.457E-01 
.449 E— 0 1 
. 365E-01 
. 1 74E-01 
.708 E-01 
. 146 E+00 
. 140 E+00 



. 1 72E-C2 
. 169E-02 
.362E-05 
.788E-C3 
. 221 E-02 
. 329E-C3 
.721E-C3 
. 434E-03 
.439E-03 
. 129E-C3 
. 2C2E-02 
. 242 E-02 
. 447E-02 
. 80 1 E-02 
.6 10 E-02 
.528 E-02 
.272E-02 
.866 E-02 
. 1C5E-01 
. 9 1 9 E-02 
. 127E-01 
. 122E-01 
•904E-02 
.796 E-02 
. 622 E-02 
.354E-03 
.821 E-02 
. 61 3E-02 
. 673 E-02 
.756 E-02 
.408 E-02 
.342 E-02 
. 209E-02 
.384 E-02 
.533E-03 
. 247 E-02 
. 126E-C3 
. 103E-02 
.680E-03 
. 180 E-02 
. 129E-C2 
. 1C3E-02 
. 254E-03 
.445E-03 
.370E-02 
.326E-03 
. 2 1 5 E-02 
. 437E-02 
. 739 E-02 
.643E-02 



.651E-04 
■ . 467 E-02 
-. 432E-03 
-.578E-02 
-. 148E-01 
-.818 E-02 
-. 904E-02 
-.550E-02 
-. 124E-02 
.791 E-02 
. 652E-02 
. 127E-01 
. 120 E-01 
. 242 E-02 

• 1 49E-0 1 
.256E-01 
. 467E-01 
.316E-01 
.272E-01 
.346E-01 
. 196E-01 
. 1 69E-0 1 
. 253E-01 
.266 E-01 
.310 E-01 
.556E-01 
.593E-01 
.401 E-01 
. 338 E-01 
.241 £-0 1 
.307 E-01 
.281 E-01 

• 467 E-0 1 
.561E-01 
. 462 E-0 1 
.332E-0 1 
. 387 E-0 1 
. 329E-0 1 
. 366 E-0 1 
.245 E-01 
.212 E-0 1 
. 170E-01 

• 1 68E-0 1 
. 964 E-02 
.21 4E-0 1 
. 484E-02 
.984 E-02 
. 171 E-01 
. 31 0E-01 
. 339E-0 1 



■ . 999E-02 
. 148E-01 
•. 167E-01 
-.191 E-01 
■ . 228E-0 1 
-.230E-01 
■ . 223E-0 1 
•. 197E-01 
-. 1 53E-01 
-.811 E-02 
■. 185E-02 
. 575E-02 
. 1 22E-0 1 
. 149E-01 
. 208E-0 1 
.286E-0 1 
.410E-01 
. 460 E-0 1 
. 478E-0 1 
.506E-01 
.47 2 E-0 1 
.426E-01 
.405E-01 
.382E-01 
.369 E-0 1 
.426E-01 
. 478E-C 1 
.456E-01 
.4 12E-0 1 
• 340E-0 1 
♦297E-0 1 
.251 E-01 
. 267E-0 1 
. 308E-0 1 
.311 E-0 1 
•273E-01 
. 258E-01 
. 228E-0 1 
.214E-01 
.1 68E-01 
. 120E-01 
.708E-02 
. 320E-02 
-. 186E-02 
• . 208 E-02 
- . 688 E-02 
-.891 E-02 
■ .793E-02 
-. 244E-02 
. 336E-02 



. 64 5E + 02 
. 645E + C2 
.645E+02 
.645E+02 
. 646E+0 2 
.646E+C2 
.646E+02 
.646E+02 
.646E+02 
•645E+02 
. 643E+0 2 
.64 1 E + 02 
. 639E+C2 
. 640E+C2 
.64 1E+02 
.640E+02 

• 640E + C2 
.643E+02 
. 643E + C2 
.643E+02 
• 647E+G2 
. 646E + C 2 
.647E+02 
.647E+02 
. 647E+C2 
. 656E+C2 
.655E+02 

♦ 661 E + 02 
.661E+C2 
. 66CE + 02 
.660E+02 
. 659E+02 
. 665E + C 2 
.665E+02 
. 667 E+02 
.669E+02 
. 67 0E + C2 
. 670E+02 
.67 0E + 02 
.671 E+02 
.670E+C2 
•669E+02 
. 668E + 02 
. 668E + C2 
.670E+02 
.674E+02 
. 674E+02 
.675E+C2 
. 677E+02 
. 675E + C2 



A-15 



TIME XNT X IT 



XI XI E 



400. 


,0 ■ 


195E+01 


402. 


,0 


. 2 6 5 E + 0 1 


404. 


,0 ■ 


-.9916+00 


406. 


,0 ■ 


- .255E+01 


408. 


,0 


•266E+01 


410. 


,0 


.126E+01 


412. 


,0 ■ 


-.413E+01 


414. 


.0 


.91 1E+00 


416. 


,0 


. 1 75E+01 


418. 


,0 


.284E+00 


420. 


.0 


. 139E+01 


422, 


.0 


.2646+01 


424. 


,0 


. 1 1 0E+C1 


426. 


,0 


. 1 66E + 0 1 


428. 


.0 ■ 


-.84 2E+00 


430. 


.0 • 


-. 159E+01 


432. 


,0 


. 1 08E + 0 1 


434. 


,0 


.103E+01 


436, 


.0 


.231 E+00 


438. 


.0 • 


-. 328E+01 


440. 


.0 


. 580E-01 


442. 


.0 


-. 136E + 00 


444 , 


.0 


-.1426+00 


446. 


.0 ■ 


- . 203E + 00 


448. 


.0 


. 141 E + C1 


450. 


.0 


-.319E+01 


452. 


.0 


-.159E-01 


454, 


.0 


- . 334E+0 1 


456. 


.0 


.133E+01 


458. 


.0 


-. 186E+00 


460. 


.0 


.591 E+00 


462, 


.0 


-.551 E+00 


464. 


.0 


- .73 9 E+00 


466, 


.0 


.1306+00 


468, 


.0 


.143E+01 


470. 


.0 


. 1 86E+00 


472. 


.0 


-.214E+01 


474. 


.0 


. 178E + 0 1 


476, 


.0 


-.490E+01 


478, 


.0 


. 1 90E+0 1 


480. 


.0 


-.4266+00 


482, 


.0 


.248 E+0 1 


484, 


.0 


-. 179E+01 


486. 


.0 


-.3836+00 


488. 


.0 


-.321 E+00 


490, 


.0 


•234E+01 


492. 


.0 


. 1 1 6E+01 


494, 


► 0 


-. 636E+C0 


496, 


.0 


. 1 89 E+0 1 


498, 


.0 


-. 139 E+0 1 


500. 


.0 


. 106E+01 



* 590E+02 
. 590E+02 
. 590E+02 
. 590E+02 
. 590E+02 
. 590E+02 
. 590E+02 
.5906+02 
.590E+02 
.590E+02 
. 590E+02 
.5905+02 
.5906+02 
. 590E+02 
. 590E+02 
.5905+02 
. 590 £+02 
.590E+02 
. 590E+02 
. 590E + 02 
.590E+02 
. 590E+02 
.590E+02 
. 590E+02 
.591 E+02 
.591E+02 
.591 E+02 
.591 E+02 
. 59 1 E+02 
.591 E+02 
. 591 E+02 
.591 E+02 
. 59 1 E+02 
.591 E+02 
.591 E+02 
.591 E+02 
. 591 E+02 
.591 E+02 
.591E+02 
.591 E+02 
.591 E+02 
. 59 1 E+02 
.59 1 E+02 
.591 E+02 
. 59 1E+02 
. 59 1 E+02 
.59 1 E+02 
.591 E+02 
.591 E+02 
.591 E+02 
.59 1 E+02 



.594E+02 
.602E+02 
. 597E+02 
. 588E+02 
•596E+02 
.599E+02 
.585E+02 
. 58 8E + 02 
. 594E + 02 
. 594E+02 
. 59 7E+02 
. 60 4E + 02 
.6046+02 
. 606E+02 
. 600E+02 
. 593E+02 
• 595E+02 
. 59 6E+02 
. 595E+02 
. 583E+02 
. 584E+02 
.584E+02 
.584E+02 
. 584E+02 
.589E+02 
.579E+02 
. 58 1 E+02 
•573E+02 
. 580E+02 
• 582E+02 
. 58 6E+02 
. 586E+02 
.585E+02 
.588E+02 
. 593E+02 
. 594E+02 
.588E+02 
. 595E+02 
. 580E+02 
. 588E + 02 
.588E+02 
. 597E+02 
.5916+02 
.590E+02 
.590E+02 
.597E+02 
.6006+02 
. 597E + 02 
. 60 2E+02 
. 596E+02 
. 598E+02 



. 403E+00 
. 1 17E+01 
.6906+00 
-. 204E+00 
. 632E+00 
. 896E+00 
-. 537E+00 
-. 200E+00 
. 34 1 E +00 
.346E+00 
. 6 83E + 00 
. 134E+01 
. 1 40E+01 
. 161E+01 
. 1 00E+01 
. 255E+00 
• 440E+00 
.581E+00 
• 455E+00 
-. 7 2 IE +00 
-. 674E + 00 
-.6776+00 
-.6656+00 
-.656E+00 
-. 142E + 00 
-. 1 10E+01 
-.926E+00 
-. 1 77E+01 
-. 1 03E+01 
-. 859E+00 
-. 463E+00 
-. 473E+00 
-. 535E+00 
-. 3 1 8E+00 
. 262E+00 
. 360E+00 
-. 275E+00 
.381E+00 
-. 1 11E+01 
-. 286E+00 
-.3036+00 
.5676+00 
-. 171E-01 
-. 8 55E-0 1 
-. 1 28E+00 
. 642E+00 
.9186+00 
. 585E+00 
. 106E+01 
. 446E+00 
. 666E+00 



A-16 



X2T 

. 127E-02 
. 127E-02 
. 127E-02 
. 1 27E-02 
. 127E-02 
. 1 27E-02 
. 1 27E-02 
. 127E-02 
. 1 27E-02 
. 1 27E-02 
. 1 27E-02 
. 127E-02 
. 1 27E-02 
. 127E-02 
.1276-02 
. 1 27E-02 
. 1 27E-02 
.1276-02 
. 1 27E-02 
. 1 27E-02 
. 127E-02 
. 127E-02 

• 1 27E-02 
. 127E-02 
. 1 27E-02 

* 1 27E-02 
.127E-02 
. 1 27E-02 
.127E-02 
. 1 27E-02 
. 1 27E-02 
. 127E-02 
. 1 27E-02 
. 1 27E-02 
. 127E-02 
. 127E-02 
. 1 27E-02 
. 1 27E-02 
. 127E-02 
. 1 27E-02 
. 1 27E-02 
. 1 27E-02 
. 127E-02 
. 127E-02 
. 1 27E-02 
. 1 27E-02 
. 1 27E-02 
. 1 27E-02 
. 1 27E-02 
. 127E-02 
. 1 27E-02 



X2 X2E 



X3 AVXNT AVX1E 



.606E-01 
. 124E+00 
. 665E-01 
230E-0 1 
. 520E-0 1 
. 697E-01 
691 fc-Cl 
-.341 E-C 1 
. I68E-C1 
. 143E-01 
. 423E-01 
. 9 64 E-C- 1 
.91 4E-01 
. 993E-01 
.321E-01 
-.419E-C1 
-. 239E-01 
-. 1 1 0E-01 
-. 232E-C1 
-. 129E+00 
-. 1 13E+00 
-. 103E + 00 
-. 906E-0 1 
- . 7 94E-0 1 
-. 230E-01 
-. 105E+00 
-. 765E-C 1 
-. 143E+00 
-.582E-01 
-. 325E-0 1 
. 107E-C1 
. 128E-01 
. 966E-02 
.31 6E-C1 
.839E-01 
.869E-01 
. 224E-0 1 
. 804E-0 1 
-. 622E-0 1 
. 188E-01 
. 1 59E-01 
. 934E-0 1 
.314E-C1 
.216E-01 
. 149E-01 
. 825E-01 
. 985E-C1 
. 573E-01 
. 929E-01 
. 258E-01 
.406E-01 



. 593 E-C1 
. 122E+00 
.653E-01 
-.242E-01 
. 515E-01 
.684E-01 
-.704E-01 
-. 354E-01 
.155E-01 
. 1 30 E-0 1 
.410E-01 
.951 E-0 1 
.902 E-0 1 
.980 E-0 1 
.308 E-0 1 
-.432 E-0 1 
-.252E-01 
-. 1 23 E-0 1 
-. 245 E-0 1 
-. 130E+00 
-.1 15E+00 
-. 104E+00 
— .919E-01 
-. 807E-0 1 
-.243E-01 
-. 107E+00 
-.777 E-0 1 
-.144 E + 00 
-.5 95 E-0 1 
-.338 E-0 1 
.947 E-02 
. 1 1 6 E-0 1 
.838 E-02 
.303 E-0 1 
. 826E-0 1 
. 857E-0 1 
.211 E-0 1 
.792 E-0 1 
-. 635E-01 
.175 E-0 1 
.146 E-0 1 
.921 E-0 1 
.301 E-0 1 
.2 03 E-0 1 
.136 E-0 1 
- .812 E-0 1 
.972 E-0 1 
.560E-01 
. 9 16 E-0 1 
. 246E-0 1 
. 393 E-0 1 



.220E-02 
. 486E-C2 
. 1 83E-02 
- . 239 E-02 
. 1 26 E-02 
.191 E-02 
-.455E-02 
- . 255 E-02 
-. 167E-04 
-. 127E-G3 
. 1 15 E-02 
. 350E-02 
. 296E-02 
. 305 E-02 
-.276E-03 
- . 360 E-02 
-. 246E-02 
-. 1 65E-C2 
-.206 E-02 
- . 667 E-02 
-. 535 E-02 
-.438E-C2 
- . 344 E-02 
-. 262 E-02 
. 174E-03 
- .358 E-02 
-. 1 94 E-02 
-. 476E-02 
-.506E-03 
.7C4E-03 
. 260 E-02 
. 246E-02 
. 209E-02 
.290E-02 
.5C0E-02 
. 469E-02 
. 134E-02 
.385E-02 
- . 296 E-02 
•975E-03 
. 754E-03 
.420E-02 
. 101 E-02 
. 473E-03 
. 126E-C3 
. 3 1 8 E-02 
. 362 E-02 
. 142E-02 
. 290E-02 
-. 395E-03 
.310E-03 



. 239E-01 
.371 E-0 1 
. 320E-0 1 
. 1 9 1 E-0 1 
. 322E-0 1 
. 382 E-0 1 
. 178E-01 
.221 E-01 
. 305 E-0 1 
.3 18 E-01 
. 383E-0 1 
.508E-01 
. 558 E-0 1 
.6 34 E-0 1 
. 59 1 E-0 1 
. 5 1 4 E-0 1 
.561 E-01 
.607 E-0 1 
.615E-01 
. 460E-0 1 
.461 E-0 1 
. 453 E-0 1 
. 444E-0 1 
. 433 E-0 1 
. 495E-0 1 
. 349E-0 1 
. 347E-0 1 
.197 E-01 
.255 E-0 1 
. 246E-01 
. 27 1 E-0 1 
. 245 E-01 
.212E-01 
.217 E-0 1 
.278E-01 
. 284E-0 1 
. 192E-01 
. 267E-0 1 
.579E-02 
. 138 E-01 
. 120 E-01 
. 223E-0 1 
. 147E-01 
. 131 E-0 1 
. 1 17 E-0 1 
.213E-01 
.260E-01 
. 23 3 E-0 1 
.308 E-01 
.251 E-0 1 
.292E-01 



. 538E-02 
. 1 12E-01 
. 1 46 E-0 1 
. 1 35E-01 
.166E-01 
. 209E-0 1 
. 182E-01 
.171 E-01 
. 1 87E-01 
.2C3E-01 
.2355-01 
.298E-01 
. 363E-0 1 
.438E-01 
. 483E-0 1 
. 493E-0 1 
.511 E-0 1 
.536E-01 
. 554E-0 1 
. 519 E-0 1 
.485 E-0 1 
.452E-01 
. 420E-01 
. 388E-01 
. 380E-0 1 
. 329 E-0 1 
. 286E-0 1 
.206E-01 
. 159E-01 
.121 E-0 1 
. 100E-01 
.789E-02 
.553E-02 
.41 3E-02 
. 524E-02 
.676E-02 
.556E-02 
.7 16E-02 
. 241 E-02 
. 1 20E-02 
- . 774 E-04 
. 230E-02 
. 222E-C2 
. 185E-02 
. 131E-02 
. 395E-02 
. 770E-02 
. 101 E-01 
. 143E-01 
.161 E-0 1 
• 1 87E-0 1 



SMCOTH 

. 677E + 02 
. 678E+C2 
. 678E+C2 
.68 1E + 02 
• 684E+02 
. 683E+C2 
.690E+C2 
.690E+02 
. 692E+02 
.691 E + C 2 
.692E+02 
.692E+02 
. 690E + 02 
.688E+C2 
•687E+C2 
.688E+02 
•688E+02 
. 688E + 02 
. 688E + C2 
.69 1E + 02 
.691E+C2 
.690E+C2 
.69CE+C2 
.689E+02 
. 690E+C2 
. 692E+0 2 
.69 1 E + 0 2 
. 692E + 02 
.691E+C2 
.69CE+C2 
. 69CE + C2 
. 690E+02 
. 690E+C2 
. 689E+C 2 
. 690E+02 
. 690E + 02 
. 692E + C2 
.693E+C2 
.700E+02 
. 70 1 E + 02 
. 70 1 E + C 2 
.7C3E+C2 
. 704E+02 
.704E+02 
.704E+02 
.706E+C2 
.705E+02 
.705E+02 
. 705E + 02 
.7C5E+02 
•705E+02 



A“»7 



001 002 003 004. O0S 



Case d) 



P MATRICES 

























• 




























r 






pii 










PM 



iaa e&i m2 m sm 



X -SCALE - 1.00E+02 UNITS'INCH. 

Y -SCALE - L00E+01 UNITS'INCH. 

CARLOS P. SIMOES MS THESIS 

OPTIMAL FILTERING OF RADAR DATA 



A-18 



000 001 002 003 001 



Case d) 



Averages and variances 
of 

measurement noise and filter error 




L -SCALE - 1.00E+02 LJNITS/INCH. 

Y -SCALE - 1.00E+00 UNITS/INCH. 

CARLOS P. SIMOES MS THESIS 

OPTIMAL FILTERING OF RADAR DATA 



A”1 9 



000 002 004; 006 005 010 



Case d) 



Absolute values 
of 

noise measurement and fitter error 




X -SCALE - L08E+02 UNITS'INCH. 

Y-SCALE - 2.00E-f00 UNITMNCh. 

CARLOS P. SIMOES MS THESIS 

OPTIMAL FILTERING OF RADAR DATA 



A “20 



010 600 500 Ti00 300 000 



Case d) 

Performance 

index 




X “SCALE - 1.88Z+0Q LJNITS/INCH. 

Y “SCALE - Q.00E+01 LJNITS/'JfiCH. 

CARLOS P. SIMOES MS TMESIS 

OPTIMAL FILTERING 'OF RADAR DATA 



A~21 



T.00 Z0O 300 100 000 



Case b) 





















































L 






pii 

222 




M 







QOO 001 002 033 004 



X -SCALE - 1.00E+02 LJNITS'im 
Y -SCALE - J.O0E-+01 UNITS' INCH. 

CARLOS P. SIMOES MS THESIS 

OPTIMAL FILTERING OF RADAR DATA 



A~22 



000 001 002 003 001 00S 




X -SCALE - L00E+02 IHITS/irm 
Y -SCALE - 1.00E+00 UNITS'INCH. 

CARLOS P. SIMOES MS THESIS 

OPTIMAL FILTERING OF RADAR DATA 



A“23 



0X0 600 S00 Ti00 300 000 



Case b) 




X -SCALE - 1.68E+02 UNITS/ INCH. 

' f -SCALE - 2.0QE408 UNIT V INCH. 

CARLOS P. SIMOES MS TMESIS 

OPTIMAL FILTERING OF RADAR DATA 



A=24 



000 002 Q&i 005 00b 010 



Case b) 




X -SCALE - 1.08E+82 IHJTS/JNCH. 

Y* SCALE - L.08E+01 UNITS'INCH. 

CARLOS P. SIMOES MS THESIS 

OPTIMAL FILTERING OF RADAR DATA 



4 



A~25 



003 001 092 003 004 00S 



Case a) 



























































All 






hi 







000 001 08Q 083 004. 



X -SCALE - 1.G8E-+0Q UNITS/INCH. 

Y-SCALE - 1.O8E401 UNITS/ INCH. 

CARLOS P. SIMOES MS THESIS 

OPTIMAL FILTERING OF RADAR DATA 



A=26 



000 001 002 003 001 




X -SCALE - LO8E+02 UNITS'INCH. 

Y -SCALE - J.00E408 UNITS/' JNCK 

CARLOS P. SIMOES MS THESIS 

OPTIMAL FILTERING OF RADAR DATA 



002 001 00S 008 010 




X -SCALE - 1.08E+02 UNITS'INCH. 

Y-SCALE - 2.e8E-f00 UNITS'INCH. 

CARLOS P. SIMOES MS THESIS 

' OPTIMAL FILTERING OF RADAR DATA 



A-28 



010 200 300 WQ 300 000 



Case a) 




X -SCALE - 1.O0E+OQ UNITS'JNCH. 

Y-SCALE - 2.00E+01 UNITS'INCH. 

CARLOS P. SIMOES MS THESIS 

OPTIMAL FILTERING OF RADAR DATA 



A“29 



S00 T-,00 ZOO 300 100 000 











} 










p 


















. 














I . 

|l V 






pii 

£22 




ni 







000 001 002 000 004 



X -SCALE - 1.00E+0Q IMTS^INCH. 

'f-SCALE - 1.00E+01 UNITS/'INCH. 

CARLOS P. SIMOES MS THESIS 

OPTIMAL FILTERING OF RADAR DATA 



A=30 



S00 Ti00 Z00 300 100 000 



Case c) 




X -SCALE - L06E+0Q UNITS/INCH. 

7-SCALE - 1.0QE+00 UHITS/INCH. 

CARLOS P. SIMOES MS THESIS 

OPTIMAL FILTERING OF RADAR DATA 



A =31 



010 §00 S00 T>00 S0O 000 



Case c ) 




X -SCALE - 1.00E402 UNITS'' INCH, 
f -SCALE - r 2.0QE40Q UNITS'' INCH. 

CARLOS P. SIMOES MS THESIS 

OPTIMAL FILTERING OF RADAR DATA 



A~32 



000 002 004 ; 006 003 010 



Case c) 




X "SCALE - 1.08E+02 UNIT S' INCH. 

V "SCALE - Q.0QE+0 1 UNJTS/JNCH. 

CARLOS P. SIMOES MS THESIS 

OPTIMAL FILTERING OF RADAR DATA 



A=33 



\ 







